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ABSTRACT 


Two  approaches  to  finite  clement  analysis  of  the  effects 
of  rock  Joints  and  fractures  on  undergroiind  unlined,  lined,  or 

r 

reinforced  tunnel  openings  are  considered.  First,  a  three-dimensional 
ubiquitous  Joint  analysis  is  described  and  applied  to  tunnel  deqj.gn. 
Although  this  approach  allows  an  estimate  of  the  weakening  effect 
of  closely  spaced  Joint  sets ,  it  is  based  on  -the  assumption  that 
the  existence  of  a  Joint  does  not  edter  the  stress  distribution 
significantly.  Since  this  may  not  edways  be  a  satisfactory 
assumption  and  because  it  may  sometimes  be  desirable  to  model 
the  action  of  individual  Joints  occurring  in  determined  locations, 
a  new  rock  Joint  representation  is  developed  and  a  Joint  stiffness 
analysis  applied  to  a  series  of  basic  problems.  The  different 
types  of  rock  Joints  and  their  relative  significance  are 

i 

discussed  and  a  classification  system  based  on  shear  strength 
of  the  Joint,  stiffness  perpendicular  to  the  Joint  and  tangential 
stiffness  of  the  Joint  is  suggested. 

Using  a  finite  element  axisymmetric  program,  states  of  stress 
in  the  rock  due  to  single  and  multiple  rock  bolt  installation  are 
examined  and  iso-stress  maps  of  the  triaxial  compression  sones  are 
shown.  States  of  stress  in  the  vicinity  of  the  rock  bolt  anchor 
and  bearing  plate  are  investigated.  The  effectiveness  of  a  single 
rock  bolt  installation  to  strengthen  rock  Joints  of  various 
orientations  is  investigated. 


i 


In  a  very  preliminary  study,  the  theory  of  a  layered 
orthographic  shell  was  coded  using  a  finite  element  plane  strain 
representation.  An  evaluation  of  the  load  resisting  capability 
of  various  tunnel  linings  incorporating  different  types  of, 
stiffeners  was  made.  Further  aneOytical  work  is  recommended. 


ii 
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ANALYSIS  OF  STRUCTURES  IN  JOINTED  ROCK 


1.  INTRODUCTION. 

This  is  a  report  of  an  investigation  performed  under  Contract 
No.  DACA  46-67“C-0016  between  the  Corps  of  Engineers  Omaha  District 
and  Dr.  Richard  E.  Goodman  as  contractor.  It  describes  results  of  re> 
search  on  analysis  techniques  performed  during  the  period  October  6,  1966 
to  June  30,  1967.  This  work  was  not  specifically  tied  to  any  project  but 
is  an  outgrowth  of  previous  analytical  work  on  the  Piledriver  Project 
under  Contract  DA-25-066-ENG-14783  completed  June  30,  1966. 

Further ,  in  continuation  of  the  work  reported  here ,  an  effort  will  be 
made  to  match  observed  structural  performance  in  the  Piledriver  test 
with  measured  material  properties.  Thus,  many  of  the  examples  and 
illustrations  have  been  developed  with  rock  and  joint  parameters  and 
loading  conditions  representative  of  the  Piledriver  test.  Frequent 
reference  is  made  to  the  final  report  of  the  previous  contract  (1) . 

In  the  previous  contract  the  principal  investigator  studied  the 
effect  of  geological  factors  on  the  behavior  of  several  experimental  tunnel 
secti(xis  of  the  Piledriver  test.  Particular  attention  was  paid  to  the  effect 
of  closely  spaced  joint  planes  on  the  strength  of  a  tunnel  under  static  load. 
The  strengthening  effect  of  rock  bolts  and  of  composite  steel  -  concrete 
liners  were  also  considered. 


1.  Goodman,  R.  E. ,  Geological  Factors  in  Design  of  Blast  Re¬ 
sistant  Tunnels  -  Piledriver  Project,  U.S.  Army  Engineer 
District  Omaha,  Technical  Report  No.  2,  Sept.  1966. 


Because  the  previous  contract  was  of  relatively  short  duration 
and  specific  direction,  a  singde  course  was  pursued  in  the  performance 
of  the  work.  It  is  the  objective  of  this  report  to  probe  the  analysis 
possibilities  more  deeply  to  develop  a  viable  technique  for  predicting 
or  back  calculating  structmral  behavior  in  jointed  rock. 

The  report  is  divided  into  three  sections.  Section  2  deals  with* 
the  analysis  of  jointed  rock.  Section  3  is  concerned  with  the  action  of 
rock  bolt  reinforcement,  and  Section  4  develops  analytical  tools  for 
evaluation  of  steel  liners. 
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2.  ANALYSIS  OF  JOINTS  AND  FRACTURES 

2.1  GENERAL 

« 

The  importance  of  weakness  surfaces  on  the  mechanics  of  rock 
is  self  evident.  The  technology  of  prediction  and  description  of  joints 
and  fractures  has,  until  recently, exceeded  the  technology  of  analysis. 

In  this  section  two  approaches  to  analysis  of  the  effect  of  joints  and 
fractures  will  be  considered.  A  three-dimensional  ubiquitous  joint 
analysis  will  be  described  and  applied  to  tunnel  design.  Then  a  quite 
different  joint  stiffness  analysis  will  be  developed  and  applied  to  a 
series  of  basic  problems.  Finally,  the  question  of  joint  classification 
will  be  discussed. 

2.2  THE  UBIQUITOUS  JOINT  ANALYSIS 

(1)  Two-Dimensional  Analysis  -  Previous  Work 

In  previous  work  on  the  Piledriver  project ,  the  importance  of 
rock  joints  on  the  probable  behavior  of  the  test  drifts  was  stated  and  a 
method  was  devised  to  make  a  behavior  prediction.  This  method  we 
have  come  to  call  the  ubiquitous  joint  analysis.  It  involved  the  calculation 
of  the  area  around  a  tunnel  within  which  the  traction  across  a  joint  of 
fixed  orientation  and  given  cohesion  and  friction  is  suHicient  to  induce 
sliding  along  or  opening  of  the  joint.  The  area  of  possible  slip  on  joints 
was  found  to  be  vastly  larger  than  the  area  of  rock  crushing  calculated 
assuming  continuous  rock. 
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The  ubiquitous  joint  analysis  assumed  (1)  that  the  orientation 
of  joints  in  a  set  can  be  measured  with  precision  and  (2)  that  a  joint  of 
a  given  set  occurs  everywhere .  hi  fact  the  analysis  may  be  teimied 
one  of  anisotropic  strength  wherein  failure  is  governed  by  two  inde¬ 
pendent  parameters  —  stress  orientation  and  stress  intensity .  An 
important  result  of  the  first  stu(fy  was  the  clear  separation  of  the 
influence  of  (1)  the  orientation  and  (2)  magnitude  of  the  principal  stresses. 

A  joint  is  stable  which  has  traction  in  the  compressive  sense 
directed  within  the  cone  of  static  friction  centered  about  its  normal. 
Therefore ,  regions  in  a  loaded  rock  body  in  which  the  joints  are  so 
oriented  relative  to  the  stress  trajectories  will  never  slip  on  the  joints 
even  if  the  intensity  of  stresses  is  large.  As  discussed  in  the  previous 
report^  the  cone  of  static  friction  centered  about  the  normal  to  a  joint  has 
a  vertex  an^e  2  <py  where  (Pj  is  the  joint  an^e  of  friction,  if  the  principal 
stress  ratio  is  0.  The  cone  containing  the  set  of  admissable  traction 
directions  across  the  joint  has  a  smaller  vertex  angle  if  the  principal 
stress  ratio  is  negative  (i.e.  one  stress  tensile)  going  to  0  as  the 
principal  stress  ratio  goes  to  -  «.  As  the  principal  stress  ratio  in¬ 
creases  positively,  the  cone  widens  until  it  contains  the  whole  joint  at  a 
principal  stress  ratio  equal  to  ctn^  (45  +  <pj/2). 

According  to  the  stress  trajectories  developed  by  a  given  loading, 
a  boimded  or  partly  bounded  slip  zone  is  defined.  For  very  large  applied 
pressures,-  or  a  very  deep  tunnel,  the  full  zone  of  slip  is  developed  for  a 
given  joint  orientation  and  strength  equation.  At  lesser  pressures,  or  les¬ 
ser  depths,  only  a  part  of  the  slip  zone  is  developed.  The  concept  of  a 
zone  of  slip  is  similar  to  Terza^i’s  concept  of  arching  in  blocky  ground. 
Beyond  a  certain  depth  corresponding  to  development  of  the  full  extent 
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of  the  "ground  arch",  the  pressures  on  tunnel  supports  cease  to  increase. 
Similarly,  beyond  a  given  pressure,  corresponding  to  full  development 
of  the  slip  zone,  the  mass  of  loosened  rock  ceases  to  eu^ent. 

A  ccxiclusion  of  the  previous  study  was  that  the  maximum  extent 
of  the  slip  zone  is  independent  of  joint  cohesion,  but  depends  simply  on 
joint  friction  and,  of  course,  on  the  ratio  of  stresses  in  the  free  field.. 
The  extent  of  the  zone  of  slip  for  any  pressure  insufficient  to  fill  out 
the  maximum  zone  depends  on  the  joint  cohesion  but  can  be  expressed 
in  dimensionless  terms. 

The  result  of  the  ubicjuitous  joint  analysis  was  a  diagram  delimit¬ 
ing  a  zone  about  the  tunnel  in  which  joints  of  a  given  set  would  tend  to 
slip.  Whether  these  zones  would  be  unstable  depends  in  fact  on  whether 
any  joints  of  the  set  actually  occur  in  the  zone  of  slip,  i.e.  on  the  actual 
joint  spacing  relative  to  the  tunnel  size.  One  can  consider  the  calculated 
slip  zone  for  a  given  joint  as  an  influence  diagram  expressing  the  region 
around  the  tunnel  in  vhich  the  joint  has  an  influence  on  stability.  The 
application  of  the  influence  diagram  demands  close  attention  to  the  geo¬ 
logical  map  which  may  indicate  the  actual  spacings  and  surface  condi¬ 
tions  of  each  set  of  joints. 

The  stability  of  the  slip  zone  also  depends  on  the  kinematic  pos¬ 
sibility  of  slip,  i.e.  on  the  existence  of  joints  of  other  sets  in  a  position 
to  open.  It  is  not  generally  theoretically  sound  to  conclude  that  a  zone  of 
overlap  of  the  influence  diagrams  of  two  existing  sets  of  joints  is  a  zone 
of  slip  on  both  sets  simultaneously.  This  would  only  be  strictly  true  if 
the  sets  were  orthogonal.  If  the  joints  are  not  orthogonal,  static  equili¬ 
brium  of  the  joint  bounded  block  demands  change  in  the  resisting  force 
mobilizing  along  each  joint  from  the  magnitude  of  the  traction  calculated 
assuming  only  one  set  to  exist  (1) .  However,  as  a  method  of  application 


1.  Goodman,  R.  E.  and  Taylor,  R,  L. ,  Proc.  Eighth  Symp.  on  Rock 
Mechanics,  Sept.  1966  (in  press). 
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it  seems  reasonable  to  direct  close  attention  to  any  region  of  the  tunnel 
showing  overlap  of  the  influence  diagrams  computed  for  two  actually  occur¬ 
ring  Joint  sets.  For  in  the  zone  of  overlap  kinematics  and  statics  offer 
the  possibility  for  movement  of  joint  blocks . 

A  method  of  tuimel  design  can  be  developed  from  the  ubiquitous 
joint  concept.  This  will  be  discussed  later. 

(2)  The  Three-Dimensional  Ubiquitous  Joint 

A  major  weakness  of  the  analysis  presented  in  the  previous  report 
is  its  limitation  to  two  dimensions.  Joint  orientations  were  expressed  in 
terms  of  their  trace  in  the  plane  normal  to  the  tuimel.  There  is  no  question 
that  this  simplification  is  artificial  for  the  actual  joint  orientations  are  three- 
dimensional  relative  to  the  tunnel  except  for  occasional  fortuitous  cases. 

The  effect  of  the  third  dimension  will  now  be  considered. 

The  ubiquitous  joint  analysis  consists  of  finding  the  normal  and 
shear  stress  on  the  joint  of  given  orientation  at  each  point  around  the 
loaded  area  or  tunnel  and  comparing  these  values  to  the  maximum  admis- 
sable  values  corresponding  to  a  straight  line  Mohr  envelope  for  the  joint 
strength: 

T  =  Ci  +  Of  tan  (p^  (2-1) 

J  J 

La  two  dimensions,  by  substituting  for  cr  and  T,  an  explicit  formula  was 
derived  defining  the  major  principal  stress  at  the  moment  of  slip,  (r^^, 
on  a  joint  oriented  at  w  relative  to  the  maximum  principal  plane  with 
principal  stress  ratio  org/or^  =  k. 
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_ J  '  J _ 

sin  ( j2  wj  -  ) 

- ^  (1  -  k)  -  (1  +  k) 

sin  (pj 


(2-2) 


At  each  point  in  a  varying  stress  field,  k  is  defined  so  the  strength 
at  each  point  can  be  determined.  An  advantage  of  this  explicit  strength 
equation  is  that  the  maximum  dimensions  of  the  zone  of  slip  can  be 
expressed  entirely  as  a  function  of  the  principal  stress  ratio  variation 
throu^out  the  rejgion. 

In  the  case  of  three  dimensions,  expressions  for  cr  and  Tmay 
be  written  in  terms  of  the  direction  cosines  of  the  normal  to  the  joint  plane 
and  the  stress  components.  Adopting  x,y,  and  z  axes  in  the  oy,  and 
cr^  directions  respectively,  the  normal  stress  is 


<7  -  0-,;  +  m2  +  n2  tTj  +  2  m  B  Ty^ 

+  2nlT^  +  21  mT^  (2-3) 

The  resultant  traction  across  the  plane  is  R  given  by 

R  —  (£  Ojjj.  +  m  Tyx  +  n  T^x)  (2  Txy  +  m  oy  +  n  T^y)^ 

+  «T^  +  mTy,,  +  li(r/ 

The  maximum  shear  T  is  thus 

T  =  \Jb^  -  <r2  (2-4) 

The  direction  cosines  for  a  given  joint  may  be  calculated  firom  the 
bearing  and  dip  of  the  joint  relative  to  the  axis  of  a  tunnel.  We  will  adopt 
the  x,z  plane  as  the  plane  of  the  tiumel  with  the  x  axis  horizontal  and  the 
z  axis  vertical.  The  y  axis  is  in  the  direction  of  the  tunnel  axis.  Then 
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the  direction  cosines  of  a  given  joint  are 
jt  =  cos  5  sin  u 

m  =  cos  5  cos  w  (2*6) 

n  =  sin  6 

where  w  is  the  bearing  of  the  normal  to  the  plane  relative  to  the  tunnel 
axis  and  6  is  the  dip  of  the  normal  to  the  plane. 

Unfortunately  substitution  of  Equations  2-3  and  2-4  in  Equation  2-1 
does  not  yield  a  tractable  equation  which  could  be  solved  to  express 
explicitly  in  terms  of  k  for  the  three-dimensional  case  in  the  manner  of 
Equation  2-2  for  the  itwo-dimensional  case.  Thus  a  different  method  of 
approach  was  followed  to  locate  the  boundary  of  the  zone  of  slip  for  any 
given  loading.  The  normal  and  shear  stresses  on  the  plane  were  calculated 
at  each  point  of  the  region.  The  strength,  S,  of  the  plane  at  each  point 
was  calculated  from  the  known  normal  stress.  The  boundary  of  the  zone 
of  slip  was  then  drawn  as  Ihe  locus  of  points  where  the  calculated  strength 
equ^s  the  shear  stress ,  T. 

There  is  a  temptation  to  regard  the  ratio  of  S  to  T  as  a  ’’factor 
of  safety”  for  each  point.  This  is  not  useful,  however,  as  S/T  of  say 
1, 1  at  a  point  in  a  body  loaded  to  say  1, 000  psi  may  remain  a  factor  of 
safety  of  1. 1  even  when  the  loading  is  increased  to  10, 000  psi.  This 
behavior  is  due  to  the  fact  that  joint  strength  depends  on  joint  orientation 
as  well  as  stress  magnitude. 

Despite  this  difficulty,  the  maximum  extent  of  the  zone  of  slip 
can  be  found  easily.  K  is  the  boundary  of  the  zone  of  slip  at  infinite  applied 
pressure.  Ih  practise,  to  delimit  the  full  zone  of  slip,  a  pressure  of 
100,000  psi  was  applied  to  the  test  tuimel  at  the  desired  free  field  principal 


stress  ratio. 
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(3)  Comparison  of  Two-Dimensional  and  Three-Dimensional  Analysis 

The  two-dimensional  ubiquitous  joint  may  be  considered  a  special 
case  of  the  three-dimensional  joint  which  lies  parallel  to  the  tunnel  axis 
(y  axis) .  To  assess  the  importance  of  considering  the  third  dimension, 
the  magnitude  of  the  slip  zone  will  be  compared  for  a  two-dimensional 
joint  of  given  slope  in  the  tunnel  cross  section  (x,z  plane)  and  the 
family  of  three-dimensional  joints  having  the  same  trace  in  the  tunnel 
cross  section  (Fig.  2-1).  This  family  can  be  generated  by  rotating 
about  thexz  trace;  i.e.  by  finding  the  set  of  direction  cosines  of  all 
planes  whose  intersection  with  the  plane  y  =  0  is  the  straight 
z  =  Sx .  We  require  intersection  of  the  plane 

Ax  +  By  +  Cz  =  0 

with  the  plane 

y  =  0 


(2-5) 

(2-6) 


obeying  the  constraint 


z  =  Sx 

Simultaneous  solution  of  these  equations  yields  the  equation  of  the 
required  planes  as 

Ax+By-—  z  =  0 


(2-7) 


(2-8) 
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The  direction  cosines  of  the  normal  to  the  family  of  planes  therefore  obey 
the  law 

f  =  Sn  (2-9) 

Since  +  m^  +  n^  =  1 

(2-10) 

Equations  2-9  and  2-10  define  the  members  of  the  family  of  planes 
passing  through  a  line  of  slope  S  in  the  tunnel  cross  section.  To  generate 
any  member  of  the  family,  one  may  fix  m  at  the  desired  value,  and  cal¬ 
culate  £  and  n.  Table  2-1  presents  values  of  £  and  n  corresponding  to 
10°  increments  in  values  of  /?  =  cos“^  m  for  joints  at  0°,  30°,  60°, 
and  90°  as  seen  in  the  tunnel  cross  section.  /3  is  defined  as  the  angle  of 
the  joint  normal  with  the  y  axis  (tunnel  axis) . 

A  joint  considered  in  the  two-dimensional  analysis  is,  in  reality, 
a  three-dimensional  joint  with  /3  =  90°,  as  for  example  joint  plane  1 
in  Figures  2.1a  and  2.1b.  To  assess  the  effect  of  ignormg  the  third 
dimension ,  the  slip  zones  on  joints  of  a  given  trace  in  the  tunnel  cross 
section  were  calculated  for  /8  in  10°  increments  from  90°  to  0°,  i.e.  on 
all  the  joint  planes  of  Table  2-1.  The  results  are  presented  in  Figures 
2,2  to  2.9. 

hi  these  figures,  the  Kirsch  solution  to  the  plane  problem  of  stresses 
about  a  tunnel  was  used.  Since  the  conditions  are  those  of  plane  strain,  the 
longitudinal  stress,  a  ,  is  at  any  point  (v)  (a  +  c  ).  hi  the  case  of  a 

y  X  z 

plane  blast  wave  travelling  across  the  tunnel  in  the  x  direction,  confine¬ 
ment  in  the  wave  front  leads  to  not  only  a  z  applied  pressure  but  a  y  applied 
pressure  as  well,  i.e.  the  wave  exerts  a  longitudinal  pressure  on  the 
tunnel.  As  a  simplifying  assumption,  we  have  assumed  that  the  longitudinal 


TAB1£  2-1.  Direction  Cosines  for  Families  of  Joints  Having  Given  Inclinations  in  the  Tunnel  Cross  Section 

(Upward  Normal  Considered  Positive ,  See  Figure  2-13b) 
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0.500  0.493  0.470  0.433  0.383  0.322  0.250  0.171  0.087  0.000 
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applied  pressure,  Py,  is  not  concentrated  aroimd  the  tunnel.  In 
Figures  2-6  to  2-9,  the  applied  longitudinal  pressure  of  Py,  as  well  as 
the  quantity  (v)  ((t^  +  cr^)  has  been  added  to  the  stress  system  in  cal¬ 
culating  the  ubiquitous  joint  slip  zones .  Li  a  two  d.unensional  analysis , 
of  course,  all  joints  parallel  the  tunnel  axis  so  longitudinal  stress  need 
not  be  considered. 

When  no  applied  longitudinal  pressure  is  considered,  Figures  2-2 
to  2-5  show  that  the  effect  of  the  third  dimension  is  (1)  to  spread  the  in¬ 
fluence  of  joints  more  evenly  around  the  tunnel  and  (2)  to  reduce  the  mag¬ 
nitude  of  the  slip  zone  particularly  at  low  pressures.  In  no  case  does 
the  influence  of  a  joint  increase  when  the  joint  is  rotated  away  from 
parallelism  with  the  tunnel  axis .  In  all  cases ,  the  slip  zone  becomes 
insignificant  when  the  angle  of  the  joint  normal  with  the  tunnel  axis  is  less 
than  about  30°.  (The  angle  of  joint  friction  assumed  in  compiling  these 
diagrams  was  31°. )  It  can  be  concluded  that  a  joint  set  which  strikes 
nearly  athwart  the  tunnel  has  little  or  no  influence  unless  it  dips  at  a 
flat  angle. 

The  90°  joint  has  its  greatest  influence  at  about  20°  of  the  wall. 

A  60°  joint  has  two  influence  zones,  a  larger  one  at  -40^  and  a  smaller 
zone  developing  at  springline  at  the  higher  pressures.  The  30°  joint  is 
important  at  +60°  and  -  60°.  The  0°  joint  is  important  at  +40°  and  -40°. 

In  the  case  of  a  vertical  blast  or  gravity  acceleration,  the  coordinate 
axes  should  be  rotated  through  90°.  Thus  a  0^  joint  is  a  horizontal  joint 
in  the  horizontal  blast  problem  but  a  vertical  joint  in  the  gravity  or  vertical 
blast  prctolem. 

Failure  to  include  an  applied  longitudinal  pressure  in  the  blast 
problem  leads  to  the  following  conclusions,  as  evidenced  by  Figures  2-6 
to  2-9.  The  joint  striking  parallel  to  the  tunnel  axis  is,  of  course. 


1.5 


2.0 


2.5 


2.0 


2.5  1.0 


r/o 


e  Diagrams 
P2/pj  =  0.43. 


im 


cohesion  in  psi. 


ORIENTATION 


ROOF 


Contour  values  are  to  be  mu 


I 


Fig.  2- 
Horizontal  a 


ORIENTATION 


ORIENTATION 


Contour  values  are  to  be  multinlied 


Fig.  2-6.  Joint  iifluence  Diagrams 
Horizontal  applied  pressure,  P2/P1  =  0.43. 
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Fig.  2-7.  Joint  Influence  Diagrams 
Horizontal  applied  pressure,  P2/P1  =  0.43. 
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affected  in  no  way  by  the  longitudinal  pressure.  All  other  joint  sets  are 
strengthened  by  the  longitudinal  applied  pressure,  i.e.  their  influence 
diagrams  are  considerably  reduced  in  size.  Further,  joints  cease  to 
have  any  weakenii^  effect  if  their  angle  with  the  tunnel  axis  is  less  than 
40°  as  opposed  to  30°  in  the  absence  of  applied  longitudinal  pressure. 

The  0®  and  30°  joints  have  their  influence  regions  shifted  away  from 
spring  line  toward  the  roof  and  floor.  The  60°  joint  has  its  slip  zone 
spread  to  the  side  wall  below  the  springline ,  while  the  90°  joint  in¬ 
fluences  the  side  wall. 

Thus,  failure  to  consider  the  third  dimension  is  conservative. 

This  same  result  has  also  been  found  in  analyzing  the  influence  of  rock 
joint  sets  on  slope  stability. 

(4)  Enlargement  of  the  Slip  Zone  -  Effect  of  Faults 

It  is  tacitly  assmned,  in  the  examples  of  the  previous  section, 
that  the  weakness  surfaces  of  the  rock  do  not  change  the  stress  distri¬ 
bution  around  the  tunnel.  This  is  believed  to  be  a  reasonable  assumpticm 
if:  (1)  the  joints  are  but  incipient  weakness  elements  revealing  them¬ 
selves  only  at  the  point  of  failure;  (2)  they  retain  residual  strength  after 
slipping;  and  (3)  the  ratio  of  minor  to  major  applied  pressures  is  not 
small,  hi  the  previous  stu(fy,  the  effect  of  stress  redistribution  due  to 
slip  on  joints  was  calculated  and  the  migration  of  the  slip  zones  during  the 
buildup  of  stress  was  plotted  for  a  number  of  cases.  It  was  found  that  if 
the  deformability  of  the  slip  zone  after  failure  was  large,  or  if  the  shear 
stress  in  the  wave  front  was  large,  the  slip  zone  augmented  considerably. 
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Should  a  major  weakness  zone  come  within  the  region  of  influence 
of  a  tunnel ,  it  remains  to  consider  whether  or  not  the  alteration  of  the 
stress  distribution  about  the  tunnel  owing  to  the  existence  of  the  fault 
could  significantly  alter  the  zone  of  slip. 

Calculation  of  the  effect  of  a  fault  on  the  stresses  around  a  tunnel 
can  be  performed  by  finite  element  analysis.  Representation  of  the  fault 
in  the  analysis  must  be  discussed.  In  the  next  section  a  joint  element 
of  a  linkage  type  is  derived  which  can  be  used  for  this  purpose,  particu¬ 
larly  if  the  fault  is  very  thin.  If  the  fault  has  ^  appreciable  thickness  of 
soft  crushed  rock  or  gouge,  it  may  be  represented  as  a  line  of  elastic 
continuum  elements  which  are  assigned  a  lower  elastic  modulus .  This 
was  the  procedure  followed  in  this  section. 

As  even  a  thick  fault  zone  is  still  tabular  in  shape ,  to  model  a 
fault  zone  with  elastic  continuum  elements,  it  is  desirable  to  loiow  the 
limiting  length  to  width  ratio  (aspect  ratio)  which  can  be  used  without 
introducing  error .  A  series  of  small  finite  element  meshes  were  made 
for  this  purpose  incorporating  a  soft  fault  seam.  The  number  of  rows 
and  the  aspect  ratio  of  the  elements  conprising  the  fault  were  varied 
systematically.  It  was  found  that  no  significant  loss  of  precision  occurred 
as  long  as  the  aspect  ratio  was  less  than  about  three  times  the  number  of 
rows  of  elements  comprising  the  fault.  To  model  a  fault,  two  rows  of 
elements  six  times  as  long  as  wide  are  about  as  precise  as  one  row  of 
elements  three  times  as  long  as  wide . 

Figure  2-lOa  presents  contours  of  maximum  shear  stress  developed 
around  a  tunnel  with  applied  pressure,  p  horizontally  and  0.43  p  verti¬ 
cally.  The  contour  values  are  expressed  as  multiples  of  p.  Figures 
2-10b,c,d,and  e  present  maximum  shear  stress  contours  for  four  cases 
in  which  a  thick  horizontal  fault  occurs  at  varying  positions  as  shown. 

The  fault  is  100  times  as  deformable  as  the  rock. 


Maximum  shear  stress  aroimd  a  circular  tunnel  with  and  without  a  fault  zone 


24 


In  the  homogeneous  case ,  the  maximum  shear  stress  occurs  at  the 
crown  and  is  about  1.25  p  (see  Fig.  2-1  Oa) .  In  Case  I  where  a  horizontal 
fault  is  centered  a/2  above  the  top  of  the  tuimel  (see  Fig.  2-1  Ob) ,  a  ma^l- 
fication  of  shear  stress  to  about  1. 8  p  occurs  and  the  center  of  TnaxiTniim 
shear  is  shifted  by  the  fault  about  10°  from  the  crown.  A  similar  but  more 
drastic  result  is  found  when  the  fault  occurs  just  above  the  crown  (see  . 
Fig.  2- 10c)  —  the  wedge  of  rock  between  the  tunnel  and  the  fault  is  under 
intolerable  shear  stresses  exceeding  2  p.  When  the  fault  intersects  the 
tunnel,  as  in  Figures  2-lOd  and  e,  no  appreciable  shear  stress  concen¬ 
tration  results. 

The  effect  of  these  stress  changes  on  the  slip  zones  of  a  30®  joint 
is  shown  in  Figures  2-11  to  2-16,  The  shear  stress  concentration  above 
the  crown  when  the  fault  is  above  the  tunnel  leads  to  an  accentuation  of  the 
slip  zone  development  in  this  region  and  the  fault  neatly  truncates  the  slip 
zone .  The  effect  of  the  fault  is  optimum  when  it  occurs  just  above  the 
tunnel.  It  is  as  if  the  occurrence  of  the  fault  above  the  tunnel  focuses 
the  applied  pressure  on  the  wall  rock.  Ih  other  than  the  optimum  positions 
of  Case  I  and  Case  n,  the  fault  changes  the  ubiquitous  joint  slip  zone 
only  in  detail . 

(5)  Application  to  Design  of  Tunnels 

The  concept  of  a  joint  slip  zone  offers  an  approach  for  estimating 
tunnel  support  forces  required  under  a  given  acceleration.  With  respect 
to  blast  resistant  design,  assume  that  the  zone  of  slip  is  developed  by  the 
passage  of  the  stress  wave  and  the  later  accelerations  hurl  the  loosened 
mass  toward  the  opening  in  the  direction  of  the  advancing  wave  front. 

With  respect  to  selection  of  tunnel  supports,  assume  development  of  the 
slip  zone  by  stress  concentrations  during  excavation  and  its  subsequent 
acceleration  downward  by  gravity.  In  each  case  if  the  assumption  is 
reasonable ,  an  approximation  of  the  net  rock  load  to  be  supported  is  the 


FIGURE  2-11 

EFFECT  OF  FAULT  ON  SLIP  ZONE 
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Efaolt  •  £rock/ioo 


FIGURE  2-13 

EFFECT  OF  FAULT  ON  SLIP  ZONE 
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integral  over  the  tunnel  wall  area  of  the  acceleration  times  the  mass  of 
rock  in  the  slip  zone  behind  each  unit  of  wall.  The  steel  sets  or  steel  or 
concrete  liner  or  rock  bolts  must  then  be  sufficiently  strong  to  reslirt 
the  net  rock  load.  Figure  2-17  illustrates  the  essential  features  of  the 
idea.  In  this  figure,  the  concrete  must  withstand 
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where  a  is  the  blast  acceleration.  The  load  is  distributed  according  to 
the  variation  of  slip  zone  depth.  The  sound  rock  between  mg  and  the  lining 
in  the  direction  of  the  blast  must  support  a  mQ  or  the  slip  zone  will  enlarge 
upward. 

Figures  2-2  to  2-9  give  the  dimensions  of  the  slip  zone  around 
circular  tunnels  for  joints  of  various  orientation  for  an  angle  of  friction 
of  31°  and  a  free  field  stress  ratio  of  0.43.  These  figures  can  be  used 
for  design  according  to  the  above  concept.  To  calculate  the  mass  and  the 
load  variation  along  the  wall  of  the  tunnel ,  the  length  of  slip  zone  behind 
each  point  of  the  wall  must  be  known,  measured  in  the  blast  propagaticm 
direction.  In  these  figures,  a  transformation  has  been  made  which  mapped 
the  circular  tunnel  wall  onto  a  straight  line,  i.e.  the  slip  zone  has  been 
plotted  on  a  Cartesian  frame  whose  y  and  x  axes  are  respectively 
0  and  r/z. 

Since  these  figures  were  developed  for  the  condition  of  maximum 
applied  pressure  horizontally,  the  direction  of  the  rock  load  is  horizontal, 
and  wo  need  to  know  the  set  of  curves  in  the  figures  describing  true 
horizontal  lines.  (In  the  case  of  a  tunnel  vmder  gravity  loads,  the 
maximum  applied  pressure  would  most  likely  be  vertical  so  the  following 
discussion  would  apply  simply  by  turning  the  figures  through  90°)  With 
reference  to  Figure  2-18 ,  a  horizontal  line  crossing  the  tunnel  wall ,  r/a  =  1 , 
at  0  =  0j^  transforms  to  a  curved  line  given  by 
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sin  0  = 


Bin  0] 
r/a 


(2-11) 


The  rock  load  behind  any  point  of  a  lining  is  proportional  to  the 
depth,  d,  of  the  failed  zone  behind  the  lining.  A  point  at  horizontal 
depth,  d,  behind  a  point  of  the  tunnel  wall  at  0  =  0j^  has  coordinates 
r  and  0  given  by 


—  =  ( —  cos  0  -  cos  01 )  (2-12) 

a  a 

Equations  2-il  and  2-12  were  used  to  construct  the  graphical  integrator 
of  Figure  2-19'.  This  figure  may  be  overlain  on  the  slip  zone  diagrams 
of  Figures  2-2  to  2-9  to  define  the  net  rock  load  and  the  load  distribution. 

An  illustrative  example  follows . 


Example  1  -  Rock  Load  in  a  Tunnel  Due  to  a  Blast 

Three  sets  of  weakness  planes  were  logged  in  a  hypothetical  tunnel 
sectiMi  as  shown  in  Figure  2-20a.  Joint  Set  1  strikes  50°  with  respect 
to  the  tunnel  axis  and  dips  42°;  it  is  spaced  at  about  5  foot  increments. 
Joint  Set  2,  spaced  every  2  feet,  stfikes  70°  with  the  tunnel  and  dips  65°. 
A  horizontal  sheared  zone,  consisting  of  numerous  closely  spaced  hori¬ 
zontal  sheared  surfaces  will  cross  the  tunnel  30°  -  45°  below  the  spring- 
line.  The  tunnel  is  circular,  20  feet  in  diameter.  It  is  to  be  subjected 
to  a  blast  with  a  horizontally  travelling  plane  wave.  The  peak  pressure 
is  10,000  psi  and  the  peak  acceleration  50  g.  The  joint  sets  are  believed 
to  have  considerable  cohesion,  ca.  1,000  psi,  whereas  the  sheared  zone 
has  only  100  psi  cohesion.  The  rock  wei^s  165  pcf. 


E  TRANSFORMATION  OF  HORIZONTAL  LINES 
DINATE  SPACE  OF  FIGURES  2-2  TO  2-9. 


001 


on  the  rig^t  to  same  scale  as  Figures 
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Table  2-2  summarizes  the  geometric  and  strength  data  for  the 
weakness  planes.  To  make  use  of  the  ubiquitous  joint  influence  dia¬ 
grams  (Figs.  2-2  to  2-9),  the  inclination,  6,  of  each  joint  in  the  tunnel 
cross  section  and  the  angle  of  the  joint  normal  with  the  tunnel  axis,  /3, 
will  be  calculated  as  follows  (Fig.  2-21): 

/3  =  cos"^  ( sin  D  sincr)  (2-lla> 

0  =  tan“^  ( tan  D  cos  a )  (2-llb) 

where  D  is  the  dip  of  the  joint  down  from  horizontal  and  positive  if  the 
y  component  is  positive,  and  a  is  the  strike  of  the  joint  relative  to  the 
tunnel  axis  and  positive  clockwise . 

0  and  for  the  three  weakness  plane  sets  calculated  by  Equation 
2-11  are  given  in  the  last  two  columns  of  Table  2-2.  It  is  apparent  imme¬ 
diately  that  Joint  Set  2  has  no  weakening  effect  on  the  tunnel  because  0 
is  less  than  40°  (31.5°) .  hi  the  region  of  occurrence  of  the  sheared  zone, 
defined  between  horizontal  lines  piercing  the  tunnel  at  0  =  -  30°  and  -  45°, 
the  overlapping  of  the  slip  zones  for  Joint  Set  1  and  the  sheared  zone  are 
as  drawn  in  Figure  2-22.  This  figure  was  obtained  as  follows: 

I 

(a)  Since  the  joint  cohesion  is  10  x  100  psi,  its  slip  zone  for  1,000  psi 
is  required.  Trace  the  portion  of  the  slip  zone  between  =  -  30°  and  -  45° 
for  Joint  Set  1  from  Figure  2-7  (30°  joint  at  /S  =  60°  with  Py  =  p^) . 

(b)  Since  the  sheared  zone  has  cohesion  100  psi,  its  slip  zone  for 
10,000  psi  is  required.  Trace  the  portion  of  the  slip  zone  between 

0j^  =  -  30°  and  -  45°  for  the  sheared  zone  from  Figure  2-6  (horizontal 
joints  at  ^  =  90°  with  Py  =  p^) . 

(c)  Trace  the  region  common  between  (a)  and  (b) . 
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TABLE  2-2,  Weakness  Planes  in  the  Rock  at  the  Ifypothetical  Tunnel 
for  Calculation  of  Rock  Load  Due  to  Blast 


Weakness 

Type 

(<A 

(c) 

psi 

strike 

Relative 
to  Tunnel 
Axis  (0) 

Dip  of 
Plane 

(D) 

Iholination  of 
Trace  in  Tunnel 
Cross  Section 
(0) 

Angle  of 
Normal  with 
Tunnel  Axis 

m 

Joint  Set  1 

5'  spacing 

31° 

1000 

+  50° 

+  42° 

+  30.0° 

59.-0° 

Joint  Set  2 

2  '  spacing 

31° 

1000 

+  70° 

+  65° 

+  36.2° 

31.5° 

Shear  Zone ,  81°  100 

occurs  in 

rock  from  30°  -  45°  below  springline 

0° 

0° 

90.0° 

TABLE  2-3. 

Weakness  Planes  in  the  Rock  at  Second  Hypothetical  Tunnel 
for  Calculation  of  Rock  Loads  on  Tunnel  Supports 

Joint  Set 

<P 

C 

a 

D 

0  90-0 

1 

31° 

0 

50 

30 

20  70 

20 

2 

31° 

0 

0 

30 

30  60 

90 

3 

31° 

0 

50 

-40 

-28  -62 

-61 
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EXAMPLES.  FIGURE  2-21  TERMINOLOGY  AND  SIGN 

CONVENTION. 
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DEPTH 


HORIZONTAL 

LINES 


NOTE:  THE  WIDTH  OF 

EACH  "RECTANGLE" 
IS  0.1 0  «  I  FOOT 


FIGURE  2-22  ROCK  LOAD  IN  EXAMPLE  I, 


note:  the  width  of 

EACH  "RECTANGLE" 
IS  O.la  »  2  FEET 


FIGURE  2-23  ROCK  LOAD  IN  EXAMPLE  2. 
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.  (d)  Note  that  the  regions  in  question  are  bounded  by  horizontal 

lines  which  transform  to  the  solid  lines  .5a  and  .  7a  in  Figure  2-19. 

{Since  ©1  is  negative,  Figure  2-19  should  be  turned  over.) 

Integrating  Figure  2-22 ,  the  net  force  per  foot  to  be  resisted  by 
the  lining  is 

F  =  yA  -  =  (165  pcf)  (0.045  X  100  ft^)  {^^) 

S  g 

=  37, 100  pounds  /  foot 

This  force  is  distributed  fairly  evenly  over  a  two  foot  section  of  the  wall 
below  the  springline .  It  could  be  sustained  by  rock  bolts ,  providing  they 
are  sufficiently  long  that  the  anchor  lies  behind  the  slip  zone. 

Example  2  -  Rock  Load  on  Tunnel  Supports 

Table  2-3  lists  data  on  the  three  closely  spaced  joint  sets  mapped 
at  the  site  of  a  circular  tunnel  40  feet  in  diameter  (see  Fig.  2-20b) .  The 
tunnel  is  1,000  feet  deep.  Joint  Set  1  with  /3  =  20°  does  not  tend  to  slip. 

The  common  region  of  slip  of  Joint  Sets  2  and  3  will  be  determined  so  that 
the  rock  load  on  tunnel  supports  due  to  gravity  can  be  calculated. 

In  this  case,  the  major  applied  pressure,  the  x  axis,  is  vertical. 

The  y  axis  remains  the  tunnel  axis  but  x  and  z  are  rotated  through  90° 
so  that  90°  -  6,  rather  than  6,  is  required  for  entering  the  slip  zone  charts. 
The  cohesion  is  10  psi.  This  means  that  the  slip  zone  will  develop  outward 
to  the  10,000  psi  contour;  thus  virtually  the  entire  zone  of  slip  is  developed. 
Stated  another  way,  the  full  height  of  the  ground  arch  develops. 

Figure  2-23  shows  the  slip  zone  common  to  Joints  2  and  3.  The 
slip  zone  for  the  -  60°  joint  (Joint  3)  is  obtained  by  turning  Figure  2-4 
upside  down.  The  slip  zone  for  Joint  2  was  obtained  by  using  Figure  2-4 
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as  presented.  Figure  2-19,  turned  90^,  was  used  to  measure  the  height 
of  the  slip  zone  at  each  point. 

The  net  rock  load  per  foot  of  tunnel  length  is 

F  =  yA  I 

=  (165  pcQ  (42  ft^/ft)  (1) 

=  6,900  poimds/foot 

The  load  is  distributed  as  shown  in  Figure  2-23. 
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2 . 3  JOINT  STIFFNESS  ANA  LYSIS 

(1)  General 

The  method  of  the  preceding  section  allows  an  estimate  of  the 
weakening  effect  of  closely  spaced  joint  sets .  It  is  predicted  on  the 
assumption  that  the  existence  of  a  joint  does  not  alter  the  stress  dis¬ 
tribution  significantly.  This  may  not  always  be  a  satisfactory  assumption. 
Further ,  it  may  sometimes  prove  desirable  to  model  the  action  of  indi¬ 
vidual  joints  occurring  in  determined  locations ,  and  the  ubiquitous  joint 
analysis  does  not  offer  this  capability.  Thus,  in  this  section,  a  new  joint 
representation  is  presented. 

(2)  The  Prototype  Joint 

A  meaningful  representation  for  joints  in  finite  element  analysis 
must  closely  model  the  actual  prototype  characteristics .  There  are  many 
types  of  joints,  and  detailed  quantitative  data  on  the  mechanical  behavior 
of  all  types  have  not  yet  been  obtained.  However,  the  following  list  of 
characteristics  can  serve  as  a  begiiming: 

(a)  Joints  are  tabular .  hi  any  two-dimensional  representation 
they  more  closely  resemble  an  irregular  line  than  a  zone  of  some  appre¬ 
ciable  thickness. 

(b)  They  have  essentially  no  resistance  to  a  net  tension  force 
directed  in  the  noiinal . 

(c)  They  offer  high  resistance  to  compression  in  the  direction  of 
the  normal  but  may  deform  somewhat  under  normal  pressure,  particularly 
if  there  are  crushable  asperities,  contiguous  altered  rock,  apertures.  Or 
compressible  filling  material. 

(d)  At  low  normal  pressures,  shearing  stresses  along  a  joint 
produce  a  tendency  for  one  block  to  ride  up  onto  the  asperities  of  the 
other.  At  high  normal  pressures,  shear  failure  along  joints  involves 
shearing  through  the  asperities . 
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(e)  As  a  first  approximation,  the  shear  strength  may  be  repre¬ 
sented  by  a  linear  Mohr  envelope,  though  a  bilinear  Mohr  envelope  would 
be  better.  Shearing  may  be  accompanied  by  dilation,  particularly  at  low 
normal  pressures, 

(f)  Small  shear  displacements  probably  occur  as  shear  stress 
builds  up  below  the  yield  shear  stress. 

The  plane  strain  continuum  element ,  whose  derivation  is  written 
fully  in  the  Piledriver  report,  cannot  satisfactorily  model  this  prototype 
behavior.  It  is  not  a  line  element  and  in  fact,  high  eccentricity  leads  to 
inaccuracy  in  stress  determination  as  discussed  in  Section  2.2  -(4). 
Instead,  the  stiffness  of  a  special  joint  element  will  be  derived.  First, 
it  will  be  useful  to  trace  the  main  steps  of  the  derivation  of  the  constant 
strain  triangle  stiffness  as  the  joint  element  will  be  designed  to  be  com¬ 
patible  with  continuum  elements  and  will  be  added  with  them  into  the 
structural  stiffiiess  matrix  (1) . 

(3)  Continuum  Element  StiEfoess  Derivation 

The  finite  element  analysis  consists  of  replacing  a  complex 
structure  by  an  assemblage  of  small  elements.  It  may  be  derived  in 
terms  of  energy ,  in  which  case  the  finite  element  approximation  may  be 
stated  as  a  decomposition  of  the  total  potential  energy  of  a  body  into 
the  sum  of  potential  energies  of  all  component  bodies.  After  calculation 
of  the  element  stiffness,  the  structural  stiffness  matrix  is  formed  at  each 
nodal  point  in  turn  by  adding  stiffaesses  of  all  elements  containing  the 
given  nodal  point . 


*A  fuller  derivation  and  complete  description  of  all  matrix  quantities  will 
be  found  in  the  previous  report ,  cited  on  Page  1 . 

1.  The  idea  of  adding  linkage  element  stiffness  to  the  total  structural 
stiffiiess  was  developed  by  Ngo  and  Scordellis,  Jour.  Amer.  Concrete 
Inst.,  V.  64,  n.  3,  March  1967. 
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To  trace  the  steps  in  formulating  the  element  sti^ess  matrix,  the 
problem  will  be  simplified  by  considering  a  one  element  mesh  with  area  A 
as  shown  in  Figure  2-24.  A  boundary  pressure,  P,  exists  along  the  side 
of  length,  L.  The  element  is  one  unit  thick. 

The  potential  energy  of  the  element,  in  the  absence  of  bo(fy  forces  is: 


J  2  ^‘"i  ^  ^ 

vol 


-I 


w^  Pj  d  A 


area 


(2-10) 


where  Wj  is  the  tensor  of  displacement  throughout  the  element,  fo  matrix 
notation 


^  J  f  (e)*  (0)  d  V  -  J  (w)*  (P)  d  A  (2-11) 

vol  area 

As  fully  described  in  the  previous  report ,  the  displacements  (w) 
may  be  expressed  in  terms  of  the  displacements  (u)  at  the  corners  (nodal 
points)  by  a  linear  interpolation  formula 

(w)  =  {(()  {(Pq)~^  (u)  (2-12) 

and  the  strains  (e)  may  also  be  written  as  a  function  of  (u) 

(e)  =  {(f')  (u)  (2-1^ 

The  stresses  (o)  are  related  to  strain  (e)  throu^  the  stress  -  strain 
matrix  (H) 

=  (H)  ((P')  (<Pq)  ^  (u) 


(0) 


(H)  (e) 


(2-14) 
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Substituting  Equations  2-12,  2-13,  and  2-14  in  the  potential 
energy  statement  2-11,  we  obtain 


<#>  =  J  f  (u)^  (H)  (^) 

yol 

-  y  (P)  dA  (2-15) 

area 

In  the  case  of  idane  strain,  none  of  the  quantities  in  the  integrands  of 
Equation  2-15  vary  throughout  the  region,  so  for  a  structure  of  unit  thick¬ 
ness,  the  integration  yields: 

<l>  =  I  A  (u)t  [  (cp')*-  (H)  (<?•)  ]  (u) 

-  L  (u)*  (P)  (2-16) 

The  potential  energy  is  minimized  by  differentiation  with  respect  to  the 
displacement  at  each  nodal  point  in  turn ,  and  equating  to  0 

N 

y  -^  =  0 

Zy  9  un 
n=l 

This  yields 

A  [  (<Po)"^^  ('P')*  (<^)  (</>(P~^]  (“)  = 

or 

A  K  u  =  F  (2-18) 

Thus  the  element  stiffness  per  unit  area  is  the  term  in  parenthesis  in 
Equation  2-17.  It  is  also  contained  in  the  parenthesis  of  2-16. 


FIGURE  2-24  STRUCTURE  CONSISTING 

OF  triangular  continuum  element. 


FIGURE  2-25  JOINT  ELEMENT- 
WIDTH  =  0;  LOCAL  COORDINATE 
SYSTEM 
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(4)  Joint  Element  Stiffnesa 


In  analogous  fashion,  beginning  from  an  energy  equation  and 
minimizing  with  respect  to  nodal  point  displacements,  we  will  derive 
the  element  stiffness  for  the  four  nodal  point  element  of  Figure  2-25. 
This  element  has  length,  L,  but  zero  width  as  the  nodal  point  pairs 
(1,4)  and  (2,3)  have  identical  coordinates.  It  is  shown  in  a  local  co¬ 
ordinate  system  with  the  x  axis  along  the  length.  The  origin  is  at  the 
center. 

The  stored  energy,  in  such  an  element  is  due  to  the  applied 
forces  per  unit  length  acting  through  the  displacements  and  must  be 
summed  through  the  element.  Thus 


iy2 

I  y  Wi  Pi  d  X 

-h/2 


In  matrix  notation 

iv'a 

^  =  \\  (w)^  (p)  d  X 

-L/2 


where 


(w)  = 


w. 


w, 


n 


top 

1 

bottom 

^8 

top 

bottom 

w„ 

=  the  relative  displacement  vector 


and 

(P)  = 


(2-19) 


(2-20) 


(2-20a) 


the  vector  of  force  per  unit  length 


(2-20b) 


4S 


The  vector  (P)  may  be  expressed  in  terms  of  a  product  of  unit 
joint  stif&ess  and  displacement 

(P)  =  (k)  (w) 

where  (k)  is  a  diagonal  matrix  expressing  the  joint  stUfness  per  unit 
length  in  the  normal  and  tangential  directions. 

0 

•Si. 

To  clarify  the  meaning  of  unit  joint  stiffness,  consider  a  direct 
shear  test  on  a  rock  joint  of  length,  L,  and  unit  width.  At  first,  a  normal 
force  is  applied  and  the  specimen  shortens  as  the  asperities  in  the  joint  . 
deform.  The  specimen  also  shortens  elastically  in  the  solid  blocks 
above  and  below  the  joint.  Subtracting  the  elastic  deformation  from  the 
total  shortening,  we  obtain  the  joint  normal  deformation,  Wjj,  which  we 
may  plot  against  the  applied  force  per  unit  length,  Figure  2-26, 

curve  1,  shows  hypothetical  data  from  such  an  experiment.  Similarly, 
as  we  apply  tangential  force ,  we  plot  the  tangential  deformation  w^  against 
the  shearing  force  per  imit  length  as  in  Figure  2-26,  curve  2.  In  the  pre- 
peak  region,  the  slopes  of  the  two  curves  give  the  imit  normal  and  unit 
tangential  stif&ess  for  the  joint  respectively. 

Substituting  Equation  2-21  in  2-20 
L/2 

>  =  I  J  (w)^  (k)  (w)  d  X  (2-22) 

-L/2 

The  displacements  (w)  may  be  expressed  in  terms  of  the  nodal  point  dis¬ 
placements  (u)  through  a  linear  interpolation  formula.  Let  u^.  and  v^  be 
displacements  in  the  tangential  and  normal  directions  respectively  at 
nodal  point  i. along  the  bottom  of  the  joint  element 


(2-21) 


(2-21a) 


Fig.  2-26.  Data  from  a  hypothetical  direct  shear  test  on  a  rock  joint. 


Specimen  Length  =  L 
Specimen  Width  =  1 
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w 


bottom 

8 


W, 


bottom 

n 


1 

2 


1-^  0 

Jj 


2  X 

1+  ^  0 
Lf 


» 


1  + 


2x 


Ui 

«2 

J 

- — 1 

_ i 

and  along  the  top  of  the  joint  element 


Ws 


top 


% 


top 


1 

2 


,  2x 


i-¥  « 

li 


» 


1  - 


2x 


[_  ju  L 

Thus  the  relative  displacement  in  the  joint  element  is 


(w)  = 


i^gtop 


-  Ws 


bottom ' 


i-A  0  -B  0  B  0  A  0 


1 

2 


-  w, 


bottom  i 


n 


0  -AO  -B  OB  0  A 


•  -J 

«3 

^3 

«4 

Mii 

^4. 

n 

0 ! 

i 

f 

“2 

^2 

j 

«3 

> 

1  — 

^3 

"4 

.^4  , 

where 


A  =  1-^ 


B  =  1  + 


2x 


Symbolically, 


(w)  =  -  (D>  (u) 


where  (D)  and  (u)  are  defined  by  Equation  2-25 . 


(2-23) 


(2-24) 


(2-25) 


(2-25a) 
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Substituting  Equation  2-25a  in  2-22  yields 


L/2 

^  =  ly  J  (u)‘(D)*(k)(D)(u)dx 

-172 

Performing  the  multiplication  indicated  by  Equation  2-26 


(2-26) 


(D)‘  (k)  (D)  = 


-A 

0 

-B 

0 

B 

0 

A 

0 


0 

-A  i 
0 
-B 
0 
B 
0 
A 


o 

1 _ 

■-A  0  -B  0  B  6  A  0  ' 

1 - 

o 

1 _ 

0  -A  0  -B  0  B  0  A 

(D)‘  (k)  (D)  = 


kgA^  0 

ABkg 

0 

-ABkg 

0 

-a\ 

0 

0  A^kjj 

0 

ABlS 

0 

-ABk^ 

0 

ABkg  0 

0 

0 

-ABkg 

0  1 

1 

0  ABkn 

0 

B^kn 

0 

-B^kn 

0 

-ABkl, 

-ABkg  0 

-B2ks 

0 

B^kg 

0 

BAkg 

0 

1  0 

0 

0 

b\ 

0 

BAkg 

-a\  0 

-ABkg 

0 

ABkg 

0 

A^kg 

®  i 

0  -A^kg 

0 

-ABkg 

0 

ABkg 

0 

a\  ! 

In  Equation  2-26  the  only  terms  that  vary  along  the  length  are  the 
three  products  of  1  -  snd  1  +  ,  i.e.  A^,  AB,  and  B^.  There  are 

thvis  three  integrals  to  be  evaluated  as  follows 


(2-27) 
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-IV2 


(2-27a)| 


(2-27b) 


<2-27c) 


Using  the  results  of  Equation  2-27,  the  expansion  of  2-26  is 


0  =  I  L  (u)*  (K)  (u) 


(2-28) 


in  which 


K 


0 

0 

0 

-2^8 

0 

0 

2kn 

Q 

Ikn 

0 
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By  analogy  to  Equation  2-16 ,  K  is  the  joint  element  stifbiess ,  per  unit 
length.  The  factor  ~  in  Equation  2-28  will  disappear  in  taking  the  variation 

A 

of  as  in  Equation  2-17. 
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The  last  step  in  the  Johit  stiffness  derivation  is  to  place  the 
element  in  a  coordinate  system  for  the  entire  structure.  Adopting  a 
gl<*al  coordinate  system  X.Y  as  drawn  in  Figure  2-27, 


X 

y 


cos 

sin 

■  X  ■ 

-sin 

cos 

Y 

- 

where  x  and  y  are  local  coordinates  in  the  tangential  and  normal  directions 
respectively. 


(5)  Computer  Code 

The  joint  element  stiffiiess  matrix  is  constructed  for  each  joint 
element  in  a  subroutine  caUed  "joint  stiff".  The  structural  stiffiiess 
matrix  for  the  entire  system  of  blocks  and  joints  is  then  assembled  by 
adding  the  appropriate  terms  of  elements  contributing  stiffiiess,  be  they 
joints  or  continuum  elements,  at  each  nodal  point  in  turn  as  described 
fully  in  the  previous  report.  After  solution  of  the  stiffiiess  equations, 
the  joint  stresses  are  calculated  from  the  known  displacements  in  a  sub¬ 
routine  called  "joint  stress".  K  the  joint  normal  stress  is  tensile  in  any 
element,  both  k_  and  k_  are  set  equal  to  zero  for  the  element  and  the 
problem  is  repeated.  Also,  the  joint  cohesion,  joint  friction,  and  resid¬ 
ual  tangential  stiffness  (see  Fig.  2-26)  are  read  in  as  data  and  the  shear 
strength  is  calculated  for  the  indicated  normal  pressure  on  each  joint. 

If  the  joint  shear  stress  exceeds  the  shear  strength,  then  kg  is  set  equal  to 

k«,  and  the  problem  is  repeated. 

®  residual 

Any  problem  can  be  restarted  where  it  was  left  off  previously  so 
that  the  number  of  cycles  can  be  increased  gr  adually  till  a  stable  state 


is  foimd. 


1 


Fig.  2-27.  Definition  of  local  coordinate  system  (x,y)  for 
joint  elements  in  terms  of  global  system  X,Y. 
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(6)  Application  to  Several  Basic  Problems 

To  gain  experience  with  the  method  and  to  assess  its  capabilities, 
several  problems  increasing  in  complexity  were  rxm.  These  include 
(a)  Patton's  problem  —  sliding  of  a  joint  with  a  tooth;  (b)  behavior  at 
joint  intersections;  and  (c)  Trollope's  problem  —  arching  of  staggered 
blocks  over  a  trapezoidal  tunnel. 

(a)  Patton's  Problem  -  Dr.  Franklin  Patton  (1)  investigated 
the  shear  strength  of  a  joint  surface  containing  a  tooth.  At  low  normal 
pressures ,  he  found  that  failure  occurs  when  one  block  rides  up  on  the 
tooth  causing  opening  of  the  flat  regions  of  the  joint.  At  high  normal 
pressures,  the  energy  required  for  one  block  to  move  up  the  asperities 
of  the  other  was  greater  than  the  energy  required  to  shear  through  the 
tooth.  While  the  shear  strength  of  the  joint  must  be  considered  purely 
frictional  when  one  block  slides  over  the  asperities  of  the  other ,  a 
cohesional  contribution  is  added  when  the  failure  is  ascribed  to  shearing 
of  the  asperities.  The  cohesion  is  a  wall  rock  property  as  will  the 
frictional  contribution  be  in  part  when  asperities  are  sheared. 

Figure  2-28a  shows  a  finite  element  model  of  the  problem.  K 
contains  16  continuum  elements  and  4  joint  elements  (heavy  lines) .  The 
bottom  block  was  restrained  and  the  top  block  was  pulled  sideways  by 
applying  nodal  point  forces  and  boundary  conditions  as  shown.  After 
three  cycles,  the  displacements  converged  to  the  solution  shown  in 
Figure  2^28b.  Only  the  stoss  side  of  the  tooth  remained  in  contact,  the 
other  joints  opening  to  the  indicated  apertures.  Shear  stress  was  high 
in  the  tooth. 


1.  Patton,  F. ,  Multiple  Modes  of  Shear  Failure  in  Rock  and  Related 
Materials,  Ph.D.  thesis.  Geology,  Univ.  of  Illinois,  1966.  (See  also 
Horn  and  Deere,  Geotechnique,  v.  Xn,  p.  319-335,  1962.) 
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(b)  Behavior  at  Joint  Intersections  -  A  basic  attribute  of  a  blocky 
rock  system  is  the  freedom  of  blocks  to  shift  and  rotate  forming  natural 
arches  and  beams.  To  reproduce  this  behavior  in  a  complex  structure, 
an  analysis  technique  must  be  able  to  simulate  block  movements  in  a  simple 
intersection.  Figure  2-29a  shows  a  mesh  used  to  simulate  the  intersection 
of  a  through  going  and  a  staggered  joint  set .  There  are  three  elastic  . 
blocks,  each  composed  of  16  elements.  The  two  joints  are  developed  by 
8  joint  elements.  The  details  at  the  Intersection  are  shown  in  Figure  2-30. 
Three  nodal  points  with  identical  coordinates  are  used  at  the  intersection. 

Two  prc^lems  were  nm.  In  Problem  1,  the  bottom  was  placed  on 
rollers  in  the  outer  corners  (A  in  Figure  2-29a).  This  creates  an  unstable 
situation  in  which  the  top  block  tends  to  drop  down  with  rotations  about  the 
Points  A.  In  Problem  2,  the  hinges  were  moved  inward  to  B  (Fig.  2-29a) 
creating  the  tendency  for  the  opposite  sense  of  rotation . 

Figure  2-31  shows  the  details  of  displacements  in  Problem  1  along 
the  joints  and  in  the  14  contiguous  elements  of  the  blocks.  In  Figure  2-31a, 
after  one  cycle,  failure  has  initiated  along  the  bottom  with  partial  opening 
of  the  vertical  joint.  After  three  cycles,  as  shown  in  Figure  2- 31b,  crack 
opening  has  propagated  upward  in  the  vertical  joint  and  has  begun  in  the 
center  of  the  horizontal  joint  whUe  the  stresses  in  the  block  corners  are 
so  high  as  to  induce  localized  crushing  or  plastic  hinging  as  inferred  by  the 
shaded  zones.  In  the  fourth  cycle,  total  collapse  was  indicated. 

Figure  2-32  shows  the  results  for  Problem  2,  with  hinges  at  B. 

As  shown  in  2-32a,  after  one  cycle,  the  vertical  joint  had  opened  at  the 
top.  After  three  cycles  the  joint  opening  had  propagated  farther  (Fig.  2-32b) 
but  the  shear  strength  in  the  horizontal  joint  had  stabilized  the  system.  Even 
though  three  joint  elements  had  failed,  the  system  reached  stable  equilibrium. 
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only  in  Problem  2. 


FIGURE  2-30. 


INTEKSECTION  OF  JOINTS.  DETAIL. 

Initially,  the  coordinates  are  the  same  for  nodal  points 
1  and  2,  for  nodal  points  4,  5  and  8  etc. 

Element  I  is  defined  by  nodal  points  2  -  5  -  4  -  1. 

Element  n  is  defined  by  nodal  points  3  -  4  -  8  -  7. 

Element  HI  is  defined  by  nodal  points  5  -  6  -  9  -  8. 
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stable  even  though  three  joint  elements  have 
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(c)  Trollope's  Problem  -  A  tunnel  in  a  system  of  staggered 
blocks.  Trollope  (1)  has  studied  the  behavior  of  trapezoidal  tunnels 
foimed  by  removing  blocks  in  a  construction  type  sequence  from  a  regular 
stacking  of  5/8"  smooth  plastic  cubes.  In  Figure  2-33,  reproduced  from 
the  cited  work,  two  regions  are  apparent  —  a  triangular  "suspended" 
zone  above  the  opening,  and  a  stable  zone  outside.  Trollope  stated  that, 

". .  .for  a  given  trapezoidal  opening,  the  'roof  can  easily  be  fretted  away 
xmtil  a  very  stable  triangular  opening  is  established".  The  simple  hori¬ 
zontal  row  of  blocks  is  a  fimdamental  unit;  it  tends  to  fail  by  shear  along 
vertical  joints  and  by  bending  as  observed  in  the  previous  example  (see 
Fig.  2-32) .  During  pure  bending,  joint  elements  in  tensile  zones  will 
initially  fail  in  tension  rather  than  in  shear. 

Figure  2-34  is  a  finite  element  mesh  to  model  the  essential  features 
of  Trollope's  block  jointed  model.  To  meet  band  width  restrictions  so 
that  a  direct  solution  method  could  be  used  in  the  computer  available 
(IBM  7094) ,  blocks  well  outside  the  expected  suspended  zone  were 
omitted  and  a  condition  of  no  horizontal  displacement  was  set  at  the  side 
margins  as  shown.  The  model  has  been  tested  with  two  different  sets  of 
values  for  the  joint  parameters  as  listed  in  Table  2-4. 


TABLE  2-4 


Joint  Parameters 

Case  I* 

Case  n* 

K 

Moderate 

High 

ks 

Moderate 

Moderate 

^®residual 

0 

2/3 

Cohesion 

0 

Moderate 

Friction  Angle 

Moderate 

Moderate 

1.  Trollope,  D.  H. ,  Felsmechanik  und  Ingenieurgeologie ,  v.  IV/3,  1966. 
*Not  to  be  mistaken  with  the  cases  discussed  in  Trollope's  paper  (1). 
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2-35.  Trollope’s  Problem  -  Case  I.  Collapse  of  Initial  Roof 
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In  Case  I,  the  initial  roof  is  unstable  as  shovra  in  Figure  2-35. 

Blocks  in  the  roof  have  slipped  out  and  rotation  of  blocks  adjacent  to 
the  opening  is  easily  recognized  in  the  figure  (the  displacements  are 
drawn  to  scale) .  The  data  used  are  from  the  third  approximation.  B;  is 
seen  from  the  figure  that  equilibrium  has  not  yet  been  achieved.  -For 
example,  the  central  block  in  the  roof  formed  is  in  the  process  of  slipping. 
Subsequently,  however,  the  roof  formed  should  be  stable. 

To  illustrate  the  behavior  of  the  model  in  Case  n,  only  the  three 
lower  block  rows  of  the  roof  are  shown  in  Figure  2-36.  With  the  joint 
parameters  mentioned,  the  initial  roof  is  stable.  Joint  elements  are 
shown  as  double  lines  when  they  have  failed  (in  tension)  and  as  heavy 
lines  when  they  are  in  compression.  To  illustrate  the  arching  effect, 
relative  values  of  compressive  forces,  where  acting  across  the  joint 
elements,  have  been  indicated  by  arrows.  The  figure  is  drawn  after  four 
approximations.  As  convergence  has  not  yet  been  reached,  the  results 
hre  not  symmetrical.  However,  it  can  readily  be  seen  that  in  approaching 
the  roof  from  above ,  the  wei^t  of  the  overlying  rock  is  being  transmitted 
to  the  abutments  via  the  outer  blocks  in  each  row .  Tensile  failure  of 
horizontal  joint  elements  lead  to  partial  separation  of  the  rows.  Thus, 
ultimately  the  main  part  of  the  lowest  block  row  will  act  as  a  beam  cariying 
only  its  own  weight.  It  is  stabilized  through  the  shear  strength  mobilized 
in  the  upper  part  of  the  vertical  joint  elements  by  the  development  of 
lateral  thrust.  This  beam  effect  decreases  in  the  subsequent  rows  above 
the  roof. 

An.  interesting  observation  (not  shown  in  Figure  2-36)  is  the  occurrence 
of  tensile  stresses  induced  by  the  lateral  thrust  in  beam  block  elements 
above  the  roof. 
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Ih  conclusion  it  should  be  emphasized  that  we  have  not  pretended 
to  discuss  all  aspects  of  Trollope's  problem,  but  have  only  used  the  model 
to  give  an  Illustration  of  the  influence  of  the  joint  parameters  on  roof 
stability.  It  seems  that  the  method  of  analysis  can  adequately  handle 
such  joint  behavior  features  as  failing  in  tension  or  shear,  rotation  of 
blocks,  development  of  arches,  and  even,  to  a  certain  extent,  the  collapse 
pattern  of  structures  in  jointed  rock. 


2.4  DISCUSSION  -  THE  DIFFERENT  TYPES  OF  JOINTS  AND  THEIR 
RELATIVE  SIGNIFICANCE 

(1)  Factors  Influencing  Joint  Parameters  -  Joint  Classification 

The  behavior  of  a  jointed  rock  mass  depends  not  only  on  the  type 
and  intensity  of  the  applied  loads  and  constraints,  but  on  (1)  the  detailed 
properties  of  component  joints  —  stiffiiess  and  strength  and  (2)  the 
fabric  of  the  joint  system  as  a  whole  —  joint  spacing,  orientation, 
persistence,  etc.  The  finite  element  analysis  procedures  described  in 
Section  2. 3  allow  these  factors  to  be  modelled  if  the  requisite  data  of  the 
prototype  characteristics  can  be  gathered.  Much  has  been  written  on  the 
methods  for  obtaining  data  on  the  fabric  of  the  joint  system  so  this  will 
not  be  discussed  here.  The  detailed  stiffness  and  strength  properties  of 
individual  joints  will  be  discussed. 

With  the  analysis  techniques  of  Section  2.3,  three  distinct  joint 
parameters  must  be  introduced:  (1)  k^,  the  unit  stiffness  across  the  joint; 

(2)  kg,  the  unit  stiffness  along  the  joint;  and  (3)  S,  the  shear  strength  along 
the  joint  (described  by  c  and  <0 .  It  is  useful  to  characterize  a  joint  by 
means  of  k„,  kg,  and  S.  This  set  of  properties  constitutes  a  necessary 
and,  for  the  time  being,  sufficient  description  of  joints  to  allow  prediction 
of  their  "potential  behavior"  under  load. 

It  is  recognized  that  the  actual  details  of  a  joint  may  include  ineffable 
or  imponderable  attributes,  and  moreover  the  characteristics  are  irregular 
and  variable  along  the  joint  (see  Fig.  2-37) .  Further  ,  such  features  as  the 
water  content  and  de^ee  of  consolidation  of  the  filling  material,  if  present, 
may  affect  all  three  joint  parameters.  While  joint  stiffiiess  and  joint 

*It  may  ultimately  prove  desirable  to  introduce  off  diagonal  stiffness  terms 
of  the  type  -kg^  to  account  for  dilation  during  shearing  in  the  analysis  of 
Section  2.3. 
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strength  may  partially  depend  upon  factors  in  common,  however,  the 
essential  independence  of  these  quantities  must  be  recognized.  The 
oft  observed  correlation  between  the  Young's  modulus  and  the  uncon¬ 
fined  compressive  strength  for  many  rock  types  should  not  lead  to  a 
misconception  that  a  similar  correlation  exists  between  joint  stifihess 

and  joint  strength.  In  an  unfilled  joint,  the  strength  is  primarily  depen^Jent 
on  the  description  of  the  vail  rock  while  the  stiffness  d^ends  primarily,, 
on  the  description  of  the  intervening  space.  The  independence  of  these 

quantities  is  illustrated  by  Figure  2-38  in  which  idealized  direct  shear 

test  curves  give  four  extreme  associations  of  tangential  stif&iess  and 

shear  strength.  It  should  be  emphasized  that  the  mode  of  failure 

(tensile,  shear)  of  a  joint  cannot  be  deduced  from  the  stiffness  values. 

The  most  important  factors  influencing  the  joint  parameters  are 
as  follows: 

will  depend  on 

-  the  contact  area  ratio  between  the  two  joint  walls 

-  the  perpendicular  aperture  distribution  and  amplitude 

-  the  relevant  properties  of  the  filling  material 
(if  present) 

kg  will  depend  on 

-  the  roughness  of  the  joint  walls  determined  by  the 
distribution,  amplitude,  and  inclination  of  asperities 

-  the  tangential  aperture  distribution  and  amplitude 

-  the  relevant  properties  of  the  filling  material 
(if  present) 

S  will  depend  on 

-  the  friction  along  the  joint 

-  the  cohesion  due  to  interlocking 

-  the  strength  of  the  filling  material  (if  present) 


Tangential  displacement 


IDEALIZED  CURVES  ILLUSTRATING  INDEPENDENCE  OF 
STRENGTH  (S)  AND  STIFFNESS  (  REPRESENTED  BY  K^). 


Fig.  2-38. 
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The  strength  of  a  joint  referred  to  in  this  context  is  the  strength  in  shear. 
The  strength  of  joints  in  tension  must  usually  be  considered  negligible. 

It  should  be  emphasized  that,  through  the  cohesion  due  to  interlocking, 
the  joint  strength  is  related  to  the  wall  rock  strength. 

Filling  material  present  in  the  joint  may  have  a  decisive  effect 
on  all  three  parameters.  Clay  material  will  in  general  indicate  low 
values  of  kj^  and  kg.  It  will  also  indicate  low  strength  values,  except 
if  an  extensive  interlocking  necessitates  shearing  of  strong  wall  rock 
asperities  during  failure.  In  particular,  the  montmorillonite  group 
minerals  can  cause  extremely  slippery  conditions  with  negligible  shear 
strength  of  plane  joints  when  the  montmorillonitic  clay  material  is  wet, 
has  a  low  consolidation,  or  is  allowed  to  swell.  Mobilized  swelling 
pressures  can  pry  the  joint  walls  apart.  If  the  joint  walls  are  coated  with 
chlorite,  talc,  graphite,  etc. ,  the  joint  will  slip  very  easily  especially  when 
wet.  Even  if  the  strength  in  these  cases  is  very  low,  the  stiffness  may  be 
high. 

A  cementation  of  the  joint,  i.e.  by  quartz,  calcite,  orepidote,  may 
give  the  joint  properties  that  are  as  good  or  even  better  than  the  proper¬ 
ties  of  the  wall  rock.  In  a  stability  context,  these  joints  may  therefore 
safely  be  ignored.  An  important  exception  is  due  to  the  possibility  of 
calcite  fillings  dissolving  with  time  especially  when  the  calcite  filling 
is  porous  or  flaky. 

Moisture  in  a  joint  will  influence  all  three  parameters  indirectly 
throu^  the  influence  on  filling  material  properties  and  will  also  directly 
influence  the  strength  of  an  unfilled  joint  through  a  reduction  in  frictional 
strength  when  the  joint  is  wet.  Water  pressures  play  as  great  a  role  in 
jointed  rock  as  in  sods;  the  effective  stress  principle  is  considered  to 
apply.  Theoretically,  pore  water  pressure  can  be  included  in  the  stiffness 


kjj  “  is  the  stiffness  perpendicular 
to  the  joint 

kg  •-  is  the  tangential  stifhvess 
S  -  is  the  joint  shear  strength 
H  -  hi^  value 
M  -  moderate  value 
L  -  low  value 


Fig.  2-39.  Classification  of  Joints  by  Means  of  1^,  kg,  and  S 
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matrix  of  the  joint.  However,  it  seems  at  present  more  relevant  to 
regard  the  pore  water  pressure  as  an  external  pressure  to  be  applied 
to  both  joint  walls.  (This  can  easily  be  done  with  the  computer  program 
used.)  Thus,  while  the  "wetness"  of  the  joint  is  considered  a  joint  properly, 
pore  water  pressure  is  regarded  as  a  system  variable. 

The  factors  influencing  the  values  of  the  joint  parameters  kjj, 
kg,  and  S  are  very  difficult  to  quantitize  for  a  given  joint.  Even  if 
relevant  data  for  each  factor  could  be  measured,  the  compilation  to  yield 
the  parameters  requires  the  dubious  task  of  deciding  on  the  relative  in¬ 
fluence  of  each  factor.  Therefore,  the  direct  measurement  of  k^,  kg, 
and  S  is  necessary. 

For  the  purpose  of  classification,  each  parameter  may  be  rated 
high,  moderate .  or  low.  This  yields  a  classification  system  of  27  joint 
parameter  combinations.  The  classification  system  is  shown  as  a  dia¬ 
gram  in  Figure  2-39.  A  joint  with  high  normal  stiffness,  moderate 
tangential  stif&iess,  and  moderate  strength  has  been  plotted  as  an  example. 

A  high  joint  stiffiiess  is  a  value  which  yields  negligible  joint  dis¬ 
placements  when  compared  to  the  elastic  displacements  of  the  rock  blocks. 
Conversely  a  low  stiffness  is  a  value  which  leads  to  joint  displacements 
that  dwarf  the  elastic  displacements  of  the  blocks.  A  moderate  joint  stiff¬ 
ness  corresponds  to  displacements  of  joints  of  the  same  order  of  magnitude 
as  the  elastic  displacements.  Similarly,  the  joint  strength  may  be  considered 
to  be  high,  moderate,  or  low  depending  on  whether  the  joints  play  a  negligible, 
participating,  or  dominant  role  in  the  strength  of  the  rock  mass. 

In  Figure  2-40,  eight  idealized  examples  of  interlocked  joints  are 
given  to  illustrate  the  influence  of  joint  aperture  and  wall  rock  strength 
on  the  joint  parameters.  The  examples  of  Figure  2-40  demonstrate  again 
that  the  joint  stiffness  and  the  joint. strength  are  independent.  The  figure 
also  illustrates  how  the  classification  system  mentioned  may  be  used.  In 
the  examples,  only  high  and  low  values  are  applied.  For  the  wall  rock 
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strength,  a  high  value  infers  that  the  resistance  to  shearing  of  the  asperi¬ 
ties  is  high.  A  low  wall  rock  strength  indicates  in  reality  highly  altered 
wall  rock. 

Example  1  -  Ti^t  joint  or  plane  cemented  joint. 

Example  2  -  As  1,  but  with  altered  wall  rock  or  plane,  tight,  and 
slippery  joint. 

Example  3  -  Perpendicularly  tight,  but  tangentially  open  joints. 

Example  4  -  As  3,  but  with  altered  wall  rock. 

Example  5  -  Tangentially  tight,  but  perpendicularly  open  joint. 

Example  6  -  As  5,  but  with  altered  wall  rock. 

Example  7  -  Open  joint  with  compressible  filling. 

Example  8  -  As  7,  but  with  altered  wall  rock  or  plane,  open,  and 
slippery  joint. 

(2)  Values  of  the  Joint  Stif&iess 

While  the  adjectives  high,  moderate .  and  low  may  be  understood 
in  a  relative  sense  when  applied  to  joint  stiffness,  the  joint  stiffness  con¬ 
cept  is  so  new  that  no  values  are  to  be  found  in  reports  and  publications 
about  joint  properties.  The  necessary  data  exist,  however,  in  the  results 
of  direct  shear  tests,  from  around  the  world,  on  specimens  with  joints. 

One  of  the  most  sophisticated  joint  test  programs  known  to  the  author 
is  a  current  investigation  in  the  Department  of  Mining  at  Imperial  College , 
University  of  London.  A  highly  controlled,  large  direct  shear  test  machine 
has  been  built  and  9-inch  width  undisturbed  specimens  of  natural  joints 
are  being  tested*.  As  an  illustration  of  the  kinds  of  stiffness  values 
associated  with  natural  joints,  the  results  of  Test  1,  performed  in  April 
1967  on  a  joint  in  quartz  porphyry  from  Rio  Tinto,  Spaiii,  will  be  discussed. 

*We  are  indebted  to  Dr.  David  Pentz  for  providing  these  data  and  allowing 
us  to  include  them  here. 
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The  joint  tested  was  9  inches  wide  by  10  inches  long.  The  angle 
of  friction  was  34°  and  the  cohesion  was  150  psi.  The  slope  of  the  pre 
peak  portion  of  the  shear  stress  -  displacement  curve  gave  a  value  of 
kg  =  3, 780  psi/inch.  The  normal  load  -  displacement  data  have  not  yet 
been  supplied  so  1<^  is  not  yet  known  to  the  writer.  The  elastic  shear 
modulus  of  the  solid  rock  was  of  the  order  of  3  million  psi  and  the  joint* 
spacing  is  of  the  order  of  5  feet.  The  elastic  shearing  displacement  of 
a  5-foot  thick  joint  bounded  block  of  rock  is  thus  given  by  a  continuum 
unit  stiffness  of  j^ches  *  3  million  psi  =  50,000  psi/inch.  Since 
the  joint  displacements  are  of  the  order  of  13  times  the  elastic  shearing 
displacements  of  the  intervening  joint  blocks,  the  joint  may  be  character- 
-iked  as  having  low  shear  stifhiess.  It  is  classified  as  a  low  stifbiess, 
moderate  strength  joint. 

(3)  The  Relative  Response  of  Different  Classes  of  Joints 

The  wide  range  of  possible  joint  conditions  indicates  the  likelihood 
of  extremely  different  response  to  applied  load  for  joints  of  different  classi¬ 
fications.  To  explore  the  possible  range  in  behavior,  a  test  problem  has 
been  selected  —  the  stresses  and  displacements  of  a  10  foot  radius  circular 
tunnel  loaded  horizontally  with  p^  and  vertically  with  P2  =  0.43  Pj^. 

The  Kirsch  solution  results  have  been  compared  with  finite  element  solu¬ 
tions  for:  (1)  a  set  of  vertic^  joints,  (II)  a  set  of  horizontal  joints,  and 
(ID)  a  system  of  vertical  and  horizontal  sets  of  joints.  In  each  system 
joint  spacing  equals  3  feet.  For  each  of  the  three  systems,  solutions  were 
obtained  for  four  extreme  cases  in  which  different  joint  properties  were 
.  ascribed  corresponding  to  the  four  permutations  of  high  and  low  joint 
stiffness.  Case  1  had  high  k^  and  high  kg  and  Case  2  had  high  kj^  and 
low  kg.  Cases  3  and  4  had  low  k^^  with  high  values  and  low  values  of 
kg  respectively  (see  Fig.  2-41) .  Selected  results  are  presented  in  Table  2-5 
and  are  summarized  as  follows . 
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(a)  Wall  Deformations  -  The  pattern  of  wall  deformations  was 
markedly  different  for  the  tunnel  in  jointed  rock  than  for  the  tunnel  in 
continuous  rock  for  each  case  in  which  one  or  both  joint  stifiEness  para¬ 
meters  was  set  low.  The  blociqr  structure  results  in  a  discontinuous 
pattern  of. deformations,  each  block  moving  independently  of  its  * 
nei^bors.  Greater  deformations  result  when  kg  is  low  than  when  k^  . 
is  set  low. 

The  case  in  which  the  joints  are  assigned  low  normal  stiffaess  and 
hi^  tangential  stiffness.  Case  3,  yields  results  quite  unlike  the  other 
three  cases  considered.  With  the  system  of  horizontal  joints,  the  wall 
moves  out  rather  than  in  and  with  the  system  of  vertical  joints  the  roof 
moves  up  rather  than  down. 

(b)  Wall  Stresses  in  the  Joint  Blocks  (9  =  qO)  -  Since  pj^  is 
horizontal,  the  point  where  8  »  0^  is  the  wall  of  the  tunnel.  The  Kirsch 
solution  at  this  free  field  principal  stress  ratio  predicts  no  tension  any¬ 
where,  even  in  the  wall.  However,  when  the  joint  systems  were  con¬ 
sidered,  in  certain  cases  significant  tension  developed  in  the  wall. 


(c)  Roof  Stresses  in  the  Joint  Blocks  (9  =  90*^)  -  The  Kirsch 
solution  predicts  a  tangential  stress  concentration  in  the  roof  (9  =  90°) 
of  2.6  pj^.  The  stresses  were  increased  greatly  for  all  three  systems 
where  the  joints  were  assigned  high  normal  stiffness  and  low  tangential 
stiffness  (Case  2) .  With  a  single  set  of  either  horizontal  or  vertical  joints 
of  this  class ,  the  maximum  stress  concentration  in  the  wall  was  5  and  6 
respectively  while  with  both  vertical  and  horizontal  joints  of  Case  2 ,  the 
joint  blocks  suffered  a  compression  stress  of  about  12  pj^.  With  a  single 
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set  of  vertical  joints  having  low  normal  stiffness  and  hi^  shear  stiffness 
(Case  3)  the  roof  was,  conversely,  completely  destressed.  With  hori¬ 
zontal  Case  3  Joints,  however,  just  the  opposite  result  was  found,  l.e. 
very  high  roof  compression  (6  p^^) . 

In  summary  this  example  shows  the  very  great  difference  in  the 
behavior  of  a  tunnel  under  given  loading  conditions  when  the  orientation 
and  properties  of  the  joints  are  varied.  Stiff,  strong  joints  do  not  alter 
the  stress  distribution  around  a  tunnel;  all  other  types  of  joints  do. 

Direct  measurement  of  joint  quantities  appears  to  be  important,  perhaps 
even  more  so  than  direct  measurement  of  the  rock  block  characteristics . 
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3.  ANALYSIS  OF  THE  ACTION  OF  ROCK  BOLTS* 

3.1  GENERAL 

In  this  section,  the  states  of  stress  due  to  single  and  multiple 
rock  bolt  installations  are  examined,  hi  most  of  the  cases ,  idealized 
rock  bolt  representations  are  used,  these  being  point  load  representations 
of  the  anchor  and  bearing  plate.  The  probable  conditions  in  the  imme¬ 
diate  vicinity  of  the  bearing  plate  and  anchor  are  investigated  in  an  attempt 
to  determine  the  error  introduced  by  the  point  load  assumption.  The 
beneficial  or  detrimental  effects  of  a  single  rock  bolt  Installed  in  a 
jointed  medium  are  investigated  and  zones  of  various  types  of  triaxial 
stress  states  are  delineated  both  for  single  and  multiple  bolt  installations. 
The  material  in  which  the  rock  bolts  are  installed  is  considered  to  be 
elastic,  isotropic,  and  continuous.  Elastic  constants,  where  required, 
are  chosen  to  be  representative  of  a  typical  rock  type  in  which  rock  bolts 
would  be  installed. 

3.2  THE  ROCK  BOLT  ANCHOR 

The  problem  of  evaluating  the  stress  distribution  in  the  vicinity 
of  a  rock  bolt  anchor  has  been  considered  by  several  investigators. 

Visual  examination  of  actual  rock  faces  which  have  been  under  anchor 
loadings,  as  well  as  examination  of  the  anchors,  reveals  that  the  stresses 
at  the  interface  between  rock  bolt  anchor  and  wall  rock  are  high  enough 
to  cause  local  plastic  yielding  of  both  substances.  Models  using  photo¬ 
elastic  techniques  have  also  helped  define  the  probable  stress  distributions. 


*This  section.was  prepared  by  Dr.  Hans  Ewoldsen. 
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If  a  mathematical  representation  of  the  single  rock  bolt  is 
desired,  the  usual  assumption  is  to  consider  both  anchor  and  bearing 
plate.  The  point  load  representation  of  the  anchor  requires  a  continuous, 
homogeneous  mediimi.  If  continuity  is  assumed,  tensile  stresses  are 
found  to  develop  behind  the  rock  bolt  anchor  which  may  exceed  the  tensile 
strength  of  the  rock.  In  an  in- situ  situation  it  is  most  probable  that  these 
tensile  stresses  do  not  develop  because  of  adjustments  on  joints  behind 
the  anchor.  The  problem  of  tensile  stresses  together  with  the  uncer¬ 
tainties  of  using  a  point  load  assiunption  indicated  that  further  study  of 
the  rock  bolt  ^chor,  as  reported  here,  would  be  desirable. 

In  order  to  model  the  presence  of  joints  in  the  physical  system ,  a 
finite  element  type  of  solution  is  used  with  which  it  is  possible  to  create 
zones  weak  in  tension.  Since  the  rock  bolt  anchor  problem  is  an 
axisymmetric  one,  the  axisymmetric  finite  element  program  was  used. 
Restrictions  inherent  in  this  approach  are  only  those  which  require  sym¬ 
metry  about  the  axis  of  the  rock  bolt.  The  total  10,000  poimd  bolt  load 
was  sustained  by  a  shearing  force  uniformly  distributed  over  die  inside 
of  the  bore  hole.  The  anchor  was  3  inches  long  and  2  inches  in  diameter. 
Assuming  a  friction  coefficient  of  0.8,  a  radial  load  of  12,500  pounds 
was  also  distributed  over  the  rock  wall. 

hi  the  first  solution,  there  were  no  weak  joints  in  the  system; 
the  stress  distributions  obtained  are  seen  in  Figures  3-la,  3-lb,  and  3-lc. 
As  can  be  seen,  tensile  stresses  develop  behind  the  anchor.  In  the  second 
solution,  a  joint  was  simulated,  occuring  immediately  behind  the  anchor 
and  perpendicular  to  the  axis  of  the  bolt.  The  resulting  stress  distri¬ 
butions  are  seen  in  Figures  3-ld,  3-le,  and  3-lf.  Li  this  case  in  the 
immediate  vicinity  of  the  anchor,  the  inclusion  of  the  joint  results  in 
stresses  being  concentrated  between  the  anchor  and  bearing  plate  with  little 
stress  being  transmitted  to  the  area  behind  the  bolt.  A  comparison  of  the 
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stress  distributions  of  the  two  cases  for  any  single  component  (radial, 
axial ,  or  tangential)  reveals  the  interesting  fact  that  the  stress  distri¬ 
butions  are  not  appreciable  different  at  points  one  or  two  anchor  dimen¬ 
sions  away  from  the  anchor.  This  would  indicate  that  it  would  make 
little  difference  in  the  zone  of  interest  lying  between  the  anchor  and 
bearing  plate  whether  a  joint  existed  behind  the  anchor  or  the  rock  was 

m 

continuous.  Thus,  the  possibility  of  tension  joints  being  created  behind 
a  line  of  rock  bolt  anchors  does  not  appear  to  be  of  great  importance  as 
the  resulting  stress  distributions  will  be  almost  i^ntical.  * 

Comparison  of  Anchor  and  Point  Load 

A  comparison  was  made  between  the  finite  anchor  and  point  load 
representations  of  the  anchor.  The  volume  of  rock  considered  in  this 
comparison  is  limited  in  extent,  thus  stress  equalities  which  would  be 
expected  at  large  radial  distances  were  not  found.  Nevertheless,  a  trend 
toward  equality  was  indicated  by  the  comparison. 

The  axial  stress  distributions  in  the  immediate  vicinity  of  the 
rock  bolt  anchor  are  shown  in  Figure  3-2  for  both  point  load  and  expansion 
anchor  representations.  The  point  load  is  applied  at  a  point  corresponding 
to  the  center  of  the  expansion  anchor.  At  points  near  the  anchor,  as  would 
be  expected,  there  is  great  difference  in  stress  distributions.  These  arise 
from  the  greatly  differing  boundary  conditions  of  the  two  solutions.  As 
radial  distance  increase,  stress  contours  of  equal  magnitude  begin  to  coin¬ 
cide  .  Thou^  the  area  included  is  not  broad  enough  to  include  areas  where 
stress  contours  from  both  solutions  would  coincide,  St.  Venant’s  principle 
and  the  trend  shown  in  the  present  study  indicate  that  coincidence  would 
occur  within  several  feet .  Thus 

*  A  zone  of  weakness  existinfr  behind  an  installation  where  closely 
spaced  rock  bolts  of  uniform  short  length  are  used  may  have 
adverse  effects.  In  this  case,  bolts  of  varying  lengths  should  be 
used. 


BOLT  LOAD  10,000  POUNDS 
TENSION  POSITIVE 


FIGURE  3-2  POINT-ANCHOR  COMPARISON. 
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in  conclusion,  the  simplifying  assvunption  of  point  loadings  is  sufficient 
provided  the  area  of  interest  is  removed  from  the  anchor.  In  the 
immediate  vicinity  of  the  anchor,  a  point  load  assumption  is  no  longer 
valid. 


3.3  ROCK  BOLT  BEARING  PLATE 


Li  most  mathematical  analyses  of  rock  bolt  stress  distributions, 
using  either  analytic  or  finite  element  solutions,  the  usual  assumption 
has  been  to  consider  the  anchor  and  bearing  plate  loadings  on  the  rock 
as  point  loads.  The  previous  section  considered  the  error  introduced 
by  replacing  the  actual  anchor  load  distribution  by  an  interior  point 
load.  This  section  considers  the  error  of  replacing  the  actual  load 
distribution  under  the  bearing  plate  by  a  point  load  at  the  surface. 

The  stress  variation  between  a  surface  point  loading  and  a  surface 
square  plate  loading  was  examined.  The  dimensions  of  the  square  plate 
(8"  X  8")  were  chosen  to  correspond  to  the  usual  bearing  plate  dimensions. 
Stresses  in  either  solution  are  linearly  dependent  upon  loading,  thus  a 
nominal  total  load  of  10,000  pounds  was  used  for  both  solutions. 

The  point  loading  solution  is  the  familiar  Boussinesq  surface 
load  on  the  half  space.  The  solution  for  a  surface  normal  load  over  a 
rectangular  domain  (uniform  load  intensity)  is  obtained  through  inte¬ 
gration  of  the  Boussinesq  solution. 

The  stresses  were  obtained  on  a  plane  containing  the  point  of  load 


application.  The  stress  variation  between  the  two  solutions  was  examined 
using  the  vertical  stress  component  0-22;  as  this  is  the  component  of 


greatest  magnitude.  The  stress  difference ,  ^22^^^^^  -  o^zzpoint’ 
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shown  in  Figure  3-3.  As  can  be  seen,  there  are  two  domains,  one 
is  which  the  stresses  of  the  point  loading  are  higher  than  the  plate,  and 
one  in  which  they  are  lower.  At  points  more  than  1-1/2  feet  from  the  point 
of  load  application,  the  difference  between  the  two  solutions  becomes 
less  than  10%.  For  most  purposes,  a  point  load  assumption  should 
result  in  sufficient  accuracy. 

3 . 4  STRESS  DISTRIBUTION  FROM  A  SINGLE  ROCK  BOLT 

The  mathematical  approximation  of  a  rock  bolt  may  be  represented 
by  the  superposition  of  an  interior  point  load  and  a  surface  rectangular 
loading.  It  was  shown  in  Section  3.3  that  for  normal  rock  bolt  dimensions, 
there  is  little  difference  between  the  plate  and  point  loadings;  therefore  in 
the  following,  the  rock  bolt  is  represented  by  the  superposition  of  normal 
surface  and  interior  point  loads.  A  nominal  load  of  10,000  pounds  was 
used.  Stresses  for  other  loadings  may  be  found  through  mxiltiplication  of 
the  10,000  pound  stresses  by  a  linear  factor. 

The  radial,  tangential,  and  axial  stress  components  (referred  to 
a  cylindrical  coordinate  system  referenced  to  the  bolt)  are  shown  in  Figures 
3-4  a,  b,  and  c  for  a  rock  bolt  10  feet  in  length. 

If  a  single  rock  bolt  is  installed,  it  is  evident  that  the  primary  bene¬ 
ficial  effect  will  be  to  increase  the  normal  stress  on  planes  perpendicular  . 
to  the  bolt  axis  between  the  anchor  and  bearing  plate .  Stresses  in  directi(ms 
normal  to  the  bolt  axis  may  be  either  tensile  or  compressive  depending  on 
position.  The  rock  bolt  stresses  were  indicated  to  be  a  linear  fimction  of 
load.  Thus  an  increase  or  decrease  in  bolt  loading  will  not  alter  the  position 
of  zero  stress  contours  but  will  result  in  change  of  intensity  of  tensile  or 
compressive  stresses. 


FIG.  3-3  POINT  -  PLATE  COMPARISON 

10,000#  total  bolt  load  10  foot  bolt  length  ^ 

Dashed  contours  give  axial  stress,  psi,  under  rock  bolt 
with  8"  X  8"  bearing  plate; 

Solid  contours  give  axial  stress  difference,  psi,  between 
solutions  for  rock  bolt  with  bearing  plate,  and  for  , 
point  load  representation# 
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FIGURE  3-4. 
LOAD  l( 


Looking  at  the  individual  stress  components,  it  is  seen  that  the 
tangential  stress  is  tensile  between  anchor  and  bearing  plate,  com¬ 
pressive  behind  the  anchor,  and  compressive  in  a  surface  lobe.  Axial 
stresses  are  compressive  between  bearing  plate  and  anchor  and 
tensile  behind  the  anchor.  The  radial  stress  distribution  is  complex 
with  lobes  of  tension  and  compression  extending  away  from  both  bear-  . 
ing  plate  and  anchor.  Stress  magnitudes  of  all  components  fall  off 
rapidly  away  from  the  bearing  plate  and  anchor. 

3.5  ROCK  BOLT  TRIAXIAL  CONDITIONS 

Superposition  of  the  respective  normal  stress  components 
allows  the  various  states  of  triaxial  stress  to  be  delineated.  These 
triaxial  stress  zones  are  presented  in  Figure  3-5.  It  is  seen  from  the 
figure  that  there  are  in  general  no  large  areas  of  triaxial  tension  or 
compression.  The  most  common  triaxial  stress  state  is  that  of  two 
compressive  and  one  tensile  normal  stresses. 

Triaxial  tension  does  exist  in  an  annular  lobe  at  the  surface . 

The  triaxially  tensile  zones  of  several  bolts  may  merge  to  form  a  surface 
"pit"  in  which  it  is  impossible  to  stabilize  material  through  rock  bolting. 
These  triaxially  tensile  zones  may  be  seen  in  rock  bolted  models  using 
crushed  rock. 

Triaxial  compression  exists  in  an  annular  lobe  near  the  surface  and 
a  very  small  lobe  near  the  anchor.  The  surface  compressive  lobe  is  of 
definite  benefit  as  lobes  of  several  bolts  may  merge  to  form  a  zone  near 
the  surface  in  which  all  normal  stresses  are  compressive;  thus  joints  will 
bo  lows  likoly  to  opon. 
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FIGURE  3-5.  VARIABLE  CONDITIONS  OF  CONFINEMENT  OWING 

TO  A  SINGLE  ROCK  BOLT. 

TOTAL  LOAD  10,000  POUNDS;  BOLT  LENGTH  10  FEET;  STRESSES  INDICATED 

ARE  TENSILE. 
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A  point  to  note  is  the  fact  that  in  the  major  portion  of  the  zone 
between  anchor  and  bearing  plate,  only  the  tangential  stress  is  tensile. 

The  radial  and  axial  stresses  are  compressive. 

3.6  THE  EFFECT  OF  JOINTS 

Ideally,  if  a  rock  bolt  is  installed  in  a  tunnel  perimeter,  the 
stresses  resulting  from  the  rock  bolt  should  increase  the  stability  of  the 
rock  surroimding  the  rock  bolt.  As  the  stress  state  of  a  single  rock 
bolt  is  a  complex  one,  there  is  no  a  priori  way  in  which  to  determine 
what  the  effect  of  the  rock  bolt  will  be  on  joint  planes  at  points  removed 
from  the  bolt. 

It  is  important  to  distinguish  the  regions  in  which  the  effect  of  a 
single  rock  bolt  is  detrimental  from  those  in  which  it  is  beneficial.  For 
this  purpose  a  ubiquitous  joint  system  (see  Section  2.2)  was  considered 
in  the  region  around  the  bolt.  A  sector  of  tunnel  perimeter  encompassing 
80®  of  arc  was  considered  with  the  bolt  length  equal  to  the  tuimel  radius. 
Joint  orientations  with  respect  to  the  tunnel  axis  were  considered  at  every 
10®  in  the  range  0®  -  90®.  Joint  planes  were  restricted  to  those  parallel 
to  the  axis  of  the  tunnel  (two  dimensional  ubiquitous  joints) . 

Using  a  criterion  of  T  =  N  tan  <j},  ranges  of  friction  angle  (^necessary 
to  prevent  slip  were  determined  at  points  in  the  rock  mass.  The  series 
of  plots  obtained  is  shown  in  Figures  3-6a  to  j.  From  these  figures  it 
is  seen  that  for  sets  of  joints  at  low  angles  to  the  bolt  axis  ,  very  high 
friction  angles  are  required  in  the  immediate  vicinity  of  the  bolt  if  slip 
is  to  be  prevented.  Zones  in  which  the  required  friction  angle  is  30°  or 
less  are  confined  to  areas  more  than  20®  of  arc  away  from  the  rock  bolt. 

As  the  joint  orientation  an^e  is  increased,  the  zones  of  low  required 
friction  angle  migrate  towar^  the  bolt  ajds ,  until  at  90®  the  low  friction 
an^e  zone  is  centered  on  the  bolt  axis. 
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Required  Friction  to  Prevent  Slip  on  Joints 

Tunnel  Radius  -  10  feet 

Single  Rock  Bolt 

Bolt  Length  -  10  feet 

Slip  Criterion  -  T  =  N  tan  <j> 

Rock  Bolt  Stresses  Only  Are  Utilized  in  Deriving  the  Required  Friction  Angle 


Explantion  of  Figures  3-6  a  -  j. 
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3.7  MULTIPLE  BOLT  STRESSES 

Single  rock  bolt  installations  are  rare;  the  more  common  usage 
finds  rock  bolts  installed  in  a  regular  patterned  arrangement,  the 
maximum  distance  between  adjacent  bolts  being  on  the  order  of  one  half 
the  bolt  length.  The  stresses  resulting  from  two  patterned  arrays  were 
examined  -  the  first  pattern  a  rectangular  array  and  the  second  a  stag¬ 
gered  array.  These  arrays  are  shown  in  Figure  3.7  and  3.8.  Loads  on  all 
bolts  were  set  at  10.000  pounds,  the  length  of  bolts  at  10  feet,  and 
spacing  at  5  feet.  Stresses  were  examined  on  cross  sections  A  -  a* 
with  the  stress  distributions  shown  in  Figures  3-7  and  3-8, 

For  the  parameter  values  chosen,  which  are  generally  repre¬ 
sentative  of  rock  bolt  installations,  the  stress  distributions  obtained 
show: 

(a)  Axial  stress  relationships  are  relatively  imchanged  by  bolt 
pattern  installation.  Superposition  of  compressive  components  from  ad¬ 
jacent  bolts  increase  the  magnitude  of  the  compressive  stresses  found  in 
the  central  two-thirds  of  the  bolt  length.  At  the  spacing  chose,  there  is 
very  httle  axial  tension  at  the  surface.  Axial  tension  exists  behind  the 
plane  of  anchors,  unless  of  course  it  is  relieved  by  joint  block  movements. 

(b)  In  the  central  portion  of  the  bolt,  both  radial  and  tangential 
tensions  are  removed  through  compressive  components  from  adjacent 
bolts.  Tensile  zones  are  retained  around  the  anchor  and  bearing  plate, 
but  the  tensHe  stress  magnitudes  of  the  single  bolt  are  reduced. 

If  the  spacing  between  bolts  is  reduced  to  2.5  feet,  the  stress 
distributions  of  Figures  3-9a.  b,  and  c  are  obtained.  It  is  seen  that  the 
stress  contour  patterns  are  simUar;  the  major  difference  resulting  from 

a  decrease  in  bolt  spacing  is  an  increase  in  the  stress  intensity  in  all 
zones  considered. 


FIGURE  3-7  STRESSES  DUE  TO  A  RECTANGULAR  PATTERN  OF  ROCK  BOLTS. 
BOLT  LOAD  10,000  POUNOSj  BOLT  LENGTH  10  FEETj  SPACING  5  FEET. 
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n’GURE  3-8  STRESSES  DUE  TO  A  STAGGERED  PATTERN  OF  ROCK  BOLTS 
BOLT  LOAD  10,000  POUNDS-,  BOLT  LENGTH  10  FEET*  SPACING  S  FEET. 
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FIGURE  3-9  STRESSES  DUE  TO  A  RECTANGULAR  PATTERN  OF  ROCK  BOLTS. 
BOLT  LOAD  10,000  POUNDS-,  BOLT  LENGTH  to  FEET;  SPACING  2.5  FEET. 
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3.8  ROCK  CONFINEMENT  CONDITIONS  RESULTING  FROM  MULTIPLE 
BOLT  PATTERNS 

Superposition  of  the  normal  stress  components  resulting  from 
a  patterne(J  array  of  bolts  results  in  a  triaxial  stress  distribution  for  the 
patterned  array.  These  rock  bolt  pattern  triaxial  zones  are  of  interest 
as  they  indicate  where  the  rock  bolts  tend  to  do  the  most  in  rock  stabilt- 
zation. 

The  triaxial  zones  existing  on  section  A  -  A'  for  the  three  cases 
of  Figures  3-7,  3-8,  and  3-9,  are  shown  in  Figures  3-10  a,  b,  and  c. 

The  most  apparent  change  from  the  single  bolt  triaxial  distribution  is  the 
development  of  a  broad  zone  of  triaxial  compression  over  the  central 
portion  of  the  bolt.  This  zone  is  seen  to  be  greater  for  the  rectangular 
array  than  for  the  hexagonal  array.  If  an  arbitrary  scale  is  set  up  in 
which  tensile  stresses  are  rated  as  minus  factors ,  axial  tension  then 
most  unfavorable ,  it  is  seen  that  the  sta^Kered  array  produces  a 

"better"  surface  zone  than  the  rectangular  array  along  the  section 
investlgaijed. 

3.9  summary 

Rock  bolts  introduce  a  complex  state  of  stress  governed  by  the 
lengths  of  the  rock  bolt  and  the  loading.  For  a  fixed  rock  bolt  geometry, 
the  stress  intensity  at  any  point  varies  linearly  with  the  loading 
magnitude . 

Triaxial  stress  zones  of  single  bolts  are  seen  to  be,  in  general, 
not  triaxially  compressive.  If  patterned  bolt  arrays  are  installed,  zones 
of  triaxial  compression  become  larger  due  to  interaction  between  adjacent 
bolts . 
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If  a  given  ficiction  angle  is  assumed,  it  is  seen  that  the  preventicm 
of  slip  on  joint  planes  is  dependent  both  on  position  and  joint  orientation. 
As  a  rule  of  thumb,  high  angle  joints  are  stabilized  near  the  bolt  axis 
and  low  angle  joints  are  stabilized  far  from  the  bolt  axis. 


4.  ANALYSIS  OF  SUPPORTS  AND  LININGS  * 

4.1  GENERAL 

Linings  and  supports  for  tunnels  are  shell  type  structures, 
in  which  rotations  as  well  as  displacements  must  be  considered. 

In  the  previous  report ,  a  ring  of  continuum  elements  was  used 
to  represent  the  lining  of  a  tunnel.  To  improve  the  representation, 
in  this  section  shell  theory  will  be  developed  and  coded  for  intro¬ 
duction  in  the  computer  solution  to  the  problem  of  support  of  tunnels 
in  jointed  rock. 


*This  section  was  prepared  by  Dr.  Robert  L.  Taylor. 
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4.2  SHELL  ELEMENT  -  PLANE  STRAIN  ANALYSIS 
(1)  General 

The  shell  segment  is  treated  according  to  layered  orthotropic 
shell  theory  (1,2,3).  hi  this  way  the  effect  of  closely  spaced  hoop  or 
longitudinal  stiffeners  placed  on  one  side  of  the  shell  can  be  incorporated 
into  the  analysis .  (Also  the  concrete  layer  can  be  combined  provided 
the  thin  shell  theory  still  applies  —  t/R  < 

The  analysis  for  the  tunnel  liner  and  the  surrounding  rock  is 
modeled,  at  present,  by  a  finite  element  plane  strain  representation. 

Thus ,  the  shell  analysis  to  follow  is  similarly  restricted  to  the  plane 
strain  condition.  A  typical  shell  is  shown  in  Figure  4-la  and  the  finite 
element  representation  in  Figure  4-lb,  hi  the  analysis  used  herein,  the 
shell  is  approximated  by  a  series  of  straight  segments  joined  at  nodes  (4,5). 
Experience  has  shown  that  provided  Lj  is  much  less  than  R,  straight  ele¬ 
ments  introduce  no  appreciable  error.  (The  additional  stored  energy  can 
easily  be  computed  as  a  measure  of  this  error.)  The  shell  and  element 
behavior  is  deduced  from  the  principle  of  minimum  potential  energy.  To 
this  end,  consider  a  typical  element  as  shown  in  Figure  4-2.  The  end 
nodes  of  the  element  are  denoted  by  I  and  J;  their  positions  are  given  by 
the  Cartesian  coordinates  x,y. 


1.  Ambartsvimyan,  S.  A. ,  Theory  of  Anisotropic  Shells,  NASA  Tech. 

Transl.  TT  F-118,  May  1964. 

2.  Dong,  S.  B.,  Pister,  K.  S. ,  and  Taylor,  R.  1.,  Jour.  Aerospace 

Sciences,  Aug.  1962. 

3.  Reissner,  E.  and  Stavslqr,  Y.,  Jour.  Applied  Mechanics ,  v.  28,  1961. 

4.  Graffton,  P.  E. ,  and  Strome,  D.  R. ,  AIAA  Jour. ,  v.  1,  Oct.  1963. 

5.  Percy,  Pian,  et.al.,  AIAA  Jour.,  v.  3,  Jan.  1965. 
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(2)  Analysis  of  Element  Stiffness 

Using  Love's  first  approximation,  the  strain  energy  density  is 
given  for  the  shell  by 


V 


[^s^s  + 


(4-1) 


where  s  is  the  shell  coordinate  and  z  is  the  normal  to  the  x  y  plane. 
For  an  orthotropic  material 


‘^s  =  «11  S  ®12  ^z 

‘^z  =  ®l2"s  °22^ 


(4-2) 


Thi-ough  the  Kirchhoff-Love  assumption  the  strains  may  be  expressed 


(4-3) 


where  tc  and  e  are  the  strains  at  some  reference  surface  in  the  shell, 

“0  *0 

i.e.  t=  0,  (In  homogeneous  shells  this  is  taken  for  convenience  at  the  shell 
middle  surface.)  and  Xs  and  Xz  are  the  change  in  curvature  of  the  shell  reference 
surface. 

Using  Equations  4-2  and  4-3,  Equation  4-1  may  be  rewritten  as 


=  i  ®11  ^s“  +  2  ®12  'z  °22  ^z^] 

V 

=  I  [  Cll  <^So  +  t  Xs)^  +  2  Ci2  (esq  +  t  Xs) 

V 

^Xz)  +  °22(%  +  txz)^]  dt  da  (4-4) 


V 
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Let  the  shell  properties  be  defined  as: 


(4-5a) 

(4-5b) 


D  .  =  C  c.,  t^  d  t 


(4-5c) 


where  h  is  the  total  shell  thickness.  Then  Equation  4-4  becomes 


V  =  5  y  [All  2  Ai2  «S|,  ‘zj  +  A22  «z„ 

A 

+  2  £g^  Xs  +  2  Bj2  (esq  Xz  +  Ezq  Xg) 

+  2  B22  °11  2  ®12  Xg  Xz 

+  D22Xz^]  da  ^  (4-6) 

For  a  cylindrical  shell  (with  an  infinite  radius  of  initial  curvature  -  flate 
plate)  the  strain  displacement  temperature  measures  are  given  by 


f|  -  “11 

(4-7a) 

€  = 
zq 

-  q!22  a  T 

(4-7b) 

CQ 

II 

02  W  ' 

"  8  S2 

(4-7c) 

Xz  = 

0 

(4-7d) 

where  A  T  is  the  temperature  change  in  the  element . 
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Thus  Equation  4-6  reduces  to 


V  •  5  J  [Aii(U_5  -  an  AT)®  -  2(Ai2U_g  + 


“22  +  DnW  gs  +  2Bii(U_s-.anAT) 

W_^J  da  *  -(4-8) 


In  order  to  maintain  compatibility  between  contiguous  sheU  elements , 
it  is  necessary  to  match  nodal  displacements  and  rotations.  To  facilitate 
this  match  U  must  have  at  least  two  degrees  of  freedom,  W  must  have  at 
least  four  degrees  of  freedom.  To  this  end  we  select  the  element  dis¬ 
placements  in  the  form 


U  =  i  [  Ui(l-P)  +  Uj(l+p)] 


(4-9a) 


W 


=  |[Wj(2-3p+p^  +  Wj(2+3p-p®) 

(1  -  p  -  p^  +  P^)  +  (-1  -  p  +  +  P^]  (4-?b) 

2  ‘  « 


where  -  1  s  p  s  1  and 


p  =  -1  .  -s 


0  :£  S  :<  L. 


(4-10) 


res- 


8Wl  , 8W 

In  the  above  0j  and  0j  are  the  nod4  rotation  j  ^  and 
pectively.  The  strain  energy  may  now  be  expressed  in  terms  of  the 
assumed  displacement  fimctions  with 


®  s=  p  H  [”j  ■ 


(4-11) 


*Some  square  terms  in  A  T  are  deleted  since  they  drop  out  when  variations 
of  V  are  taken. 
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W.ss  =  W_pp(P,/ 

«  [  Wj  •  6  p  -  Wj  •  6  p 

1 

0j  Xi 

+  ^  (6  P  -  2) 


In  order  to  match  shell  displacements  with  the  rock  displacements 
(continuum  elements),  the  shell  displacements  Uj,  Wj,  etc.,  are  express 
ed  in  terms  of  Uj  and  Vj  by  the  transformation  equation  (see  Fig.  4-2) . 


Uj 

L  "l  *  L  ’'l 

(4-12a) 

li 

-  —  Ut  +  T  Vt  (etc.  for  Uj  and  Wj) 

Xj  Xi  ■** 

(4-12b) 

Equation  4-8  is 

expressed  finally  in  terms  of  Uj,  Vj,  0j,  Uj,  Vj,  and  0j 

+  1 

as  : 

V  = 

f  -  TiSi  dp 

(4-13) 

where 

- 1 

^11  = 

^11-  *'l2  =  ®11-  ''22  =  “ll' 

(4-14a) 

It 

-^11  “11  ^  ^12  “22  ^ 

(4-14b) 

Tg  = 

BiiaiiAT  +  BjgaiiAT  +  Bjg  ^22  ^  T 

and 

Si  = 

®2  =  W  33 

(4-14C) 

The  Sj.  are  related  to  the  nodal  displacements  through 

II 

(4- 15a) 

0,L 


+  ^(6p+2)]  ..(4-llb) 


where 
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where 

“  (“r  “J’  * 

,  A  BAB  6B  p  6  A  p 
"  \"  l2  ’  ■  L2  '  l2  *  I?  ’  ”  l3  ’  ’ 

3  p  -  1  6_B_b  6  a  p  3  p  +  1  \ 

L  ’  L3  ’  ‘  l3  ’  L  ' 

The  non  zero  are 

"  ^122  ~  ®153  “  ®164  "  ^215  “  ®226 
=  ^237  "  ^258  '  ^269  '  ^2,7,10  ""  ^ 


(4-15b) 


.  (4-15C) 


(4-15d) 


Thus  Equation  4-13  becomes 

V  “  I  «jGijk^i£®^knG£mnQ!m  "  “jGijk^iYk 

=  2  “j  ®in  "  “i 

where  ^  2. 

»ta=  J  ^  dp 

-1 

+  1 

-1 


*Numbering  is  for  computer  coding  applications . 

**It  is  to  be  noted  that  the  above  development  depends  on  use  of  the  com¬ 
puter  to  perform  the  indicated  operations  and  construction  is  such  that 
Gaussian  quadriture  may  be  used  for  the  integrations.  The  number  of 
operations  does  not  excessively  surpass  the  steps  necessary  to  com 
pute  the  element  stiffness  by  hand  explicity  before  programming  is  at¬ 
tempted  (see  listing)  but  greatly  simplifies  the  coding. 
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The  kjjjj  is  the  element  stiffness  and  fj  the  thermal  load  vector.  Assembly 
into  the  shell  structure  stif&iess  proceeds  in  the  usual  manner  by  the  direct 
stiffness  method. 


(3)  Shell  Properties 

The  shell  properties  for  layered  shells  are  obtained  from  Equation 
4-5  by  changing  the  Cjj  for  each  layer.  Accordingly  one  may  write 

N  +  l  N 

=  I  I  =  1  “I 


k=l  h. 


k=l 


where  N  is  the  number  of  layers  in  the  shell,  ^  are  the  properties  of 
the  ''k"  layer,  h,  is  the  distance  to  the  most  negative  surface  of  the  k  layer, 
and  h|^^  ^  is  the  distance  to  the  most  positive  surface.  Similarly 
N 


®ij  ■  [\+  1  ■ 

k=l 


N 


k-1 

For  a  single  layer,  the  above  definitions  reduce  to  those  conventionally 
used  in  shell  theory. 
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(4)  Shell  Forces 

The  forces  in  each  element  may  be  computed  from  the  values 
of  the  nodal  displacements  and  rotation  using  the  constitutive  equations 
together  with  Equations  4-7  and  4-11.  The  constitutive  equations  con¬ 
sistent  with  Love's  first  approximations  are  given  by 


Nz  = 

Mg  = 

Mz  = 


^11  ®s  ^12 

Ai2  ^22 

®11  *^s  ®12 


BllXs  +  Bi2Xz 
Bi2Xs  +  ^22X2 
^11  ^s  +  ^2X2 
^12  ^  ^22  ^z 


For  cylindrical  shells  in  plane  strain  the  most  important  shell  forces 

are  N  and  M  which  are,  with  aid  of  Equation  4-7,  given  by 
s  s 

Ns  =  +  A^2  Xg 


®11  '^s  ®12  ^11  ^S 


Now  Equation  4-11  completely  defines  each  shell  force  over  each  element 
in  terms  of  the  nodal  displacements.  In  the  computer  codes  these  are 
evaluated  only  at  the  midlength  of  each  element.  Ih  addition  the  shear  is 
evaluated  from  the  equilibrium  requirement 

9  Mg 

as 

with  the  displacement  expansion  used.  The  above  expression  reduces  to 
Qs  =  -I^llW^sss 


118 


(5)  Computer  Code 

The  above  development  has  been  coded  in  FORTRAN  IV  language. 
Preliminary  debugging  was  accomplished  for  simple  beam  type  configur¬ 
ations,  The  program  was  then  used  for  a  comparative  study  of  the  per¬ 
formance  of  stiffened  shell  linings  as  described  below. 


4.3  SHELL  LINERS  -  EVALUATION  OF  ALTERNATE  DESIGNS 

To  assess  the  effectiveness  of  tunnel  linings  In  resisting  loads, 
five  tunnel  liners  were  analyzed.  The  sections  selected  for  analysis 
are  alternate  designs  incorporating  different  types  of  stiffeners  (see 
Fig.  4-3).  Section  A  is  an  unstiffened  steel  liner.  In  Section  B  a  5/8  ^ 
inch  plate  has  been  welded  to  the  steel  to  make  7.2  inch  deep  stiffeners 
12  inches  on  centers.  In  Section  C,  similar  stiffeners  have  been  added;  ’ 
they  are  2-1/2  inches  deep  on  6  inch  centers  and  are  formed  of  1/4  inch 
plate.  All  sections  contain  equivalent  weights  of  steel. 

The  finite  element  mesh  used  in  the  analysis  is  shown  in  Figure 
4-4.  The  liner  appears  as  the  circular  line  segment  and  has  an  inner 
radius  of  42  inches.  Under  the  loading  it  is  assumed  that  a  portion  of  the 
rock  and  concrete  liner  fail  and  will  carry  no  load.  The  failed  portion 
is  then  accelerated  against  the  liner  as  in  the  example  of  Section  2. 2- (5). 

The  analysis  utilized  a  pressure  load  on  the  shell  to  simulate  the  rock 
acceleration  load.  The  same  load  has  been  used  for  each  shell  configu¬ 
ration,  hence,  the  relative  effectiveness  of  each  liner  with  respect  to  the 
others  may  be  assessed.  During  this  preliminary  investigation,  only  the 
displacement  profiles  have  been  computed.  Also,  only  elastic  shell  behavior 
is  included.  The  cases  analyzed  are  listed  in  Table  4-1  and  described  in 
Figure  4-3. 

The  displacement  plots  for  Cases  1-Al,  1-A2,  and  1-B  are  shown 
in  Figure  4-5  and  for  Cases  1-A2  and  1-C  in  Figure  4-6.  From  these 
figures  it  is  observed  that  stiffeners  will  decrease  displacements  in  the 
liner  appreciably.  The  chances  of  complete  collapse  of  the  1-Al  and 
1-A2  configurations  by  the  formation  of  plastic  hinges  in  the  region  of 
hi^  curvature  should  be  considered  in  future  analyses. 
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TABLE  4-1 


Case 

Tunnel  Liner 
No. 

Liner  Section 

Thickness 

(in.) 

Concrete 
f'c  -  (psi) 

1-Al 

1 

A 

1.000 

3500 

2-Al 

2 

A 

1.000 

9000 

1-B 

•  3 

B 

0.625 

3500 

1-A2 

5 

A 

0.375 

3500 

1-C 

7 

C 

0.375 

3500 
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FIGURE  4-5  LINER  DISPLACEMENTS.  FIGURE  4-6  LINER  DISPLACEMENTS. 

CASES*  l-AI,  l-B  AND  2-AI.  CASES  I  -  A  2  AND  l-C. 
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